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Abstract. 

We investigate the mechanism by which topological defects form in first order phase 
transitions with a charged order parameter. We show how thick superconductor 
vortices and heavy cosmic strings form by trapping of magnetic flux. In an external 
magnetic field, intermediate objects such as strips and membranes of magnetic flux and 
chains of single winding defects are produced. At non-zero temperature, a variety of 
spontaneous defects of different winding numbers arise. In cosmology, our results mean 
that the magnetic flux thermal fluctuations get trapped in a primordial multi-tension 
string network. The mechanism may also apply to the production of cosmic-like strings 
in brane collisions. In a thin type-I superconductor film, flux strips are found to be 
meta-stable while thick vortices are stable up to some critical value of the winding 
number which increases with the thickness of the film. In addition, a non-dissipative 
Josephson-like current is obtained across the strips of quantized magnetic flux. 

Keywords: Gauge theories, superconductivity, cosmic strings. 

1. Introduction 



The Kibble mechanism [T] has proved successful in explaining the formation of 
topological defects in phase transitions of field theories which experience spontaneous 
breakdown of a global symmetry. It applies to both first and second order phase 
transitions. The vacuum (minimum energy configuration) is degenerate and the system 
can fall into any vacuum state in different spatial locations. In first order phase 
transitions, where bubbles nucleate spontaneously, it has been tested theoretically 
in numerical simulations and analytical calculations [21 El H]. Experimentally it has 
been verified in nematic liquid crystals through the spatial distribution of disclinations 

mum- 

In second order phase transitions Kibble's argument was completed by Zurek [8] 
on what is known as Kibble-Zurek mechanism (KZ). It applies to the generation of 
topological defects as the system in question cools down in a finite time. When the 
transition takes place, non-equilibrium phenomena make the coherence length finite. 
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This allows the system to choose different vacuum states in different locations. The 
mechanism has been tested experimentally in the prediction of the vorticity pattern of 
superfluid ^He ^[10]. 

In cosmology, the analogues of disclination lines and superfluid vortices are global 
cosmic strings. They may have been generated in the early stages of the universe 
[TT| [12], [13]. If they exist. Kibble mechanism is expected to explain their origin. 

The order parameters which describe the phases of liquid crystals, ^He and field 
theories with global cosmic strings are neutral. In the contrary, the order parameter 
which describes the symmetry of a gauge theory carries a charge and is coupled to a 
gauge field. KZ mechanism is based on the dynamics of a neutral order parameter, and 
therefore it is insufficient for describing generally either the formation of local cosmic 
strings in cosmological phase transitions of gauge theories or the origin of vortex lines 
in superconductors ^14] . 

It is easy to see this in a simple condensed matter example. When a type-II 
superconductor is slowly cooled from the normal to the superconducting phase in an 
external magnetic field, an Abrikosov lattice consisting of vortices of the same sign is 
formed. This cannot be explained by Kibble-Zurek mechanism because it predicts the 
same number of vortices and anti-vortices. In the same way. Kibble mechanism is not 
sufficient to explain the formation of intermediate structures such as strips of magnetic 
flux in type-I superconductors under an external magnetic field [151 HB] • 

For cosmology and particle physics prospects, gauge theories are more relevant 
than global theories. Grand Unified Theories (GUTs) are gauge theories which predict 
generically the spontaneous formation of local cosmic strings, monopoles and other 
topological defects. In particular, cosmic strings are tubes of quantized magnetic flux. 
The value of the tension of the strings increases with their flux. The phenomenological 
effects of these strings are mainly gravitational and depend quantitatively on their 
tension. Different cosmological observations generally attribute different upper bounds 
to the tension of cosmic strings [17]. These are average values which could actually 
correspond to a distribution of windings as a result of the evolution of a primordial 
multi-tension string network [18]. Kibble mechanism, on the contrary, gives rise to a 
unique tension because all strings have unit winding. 

In a more theoretical context, recent findings suggest that cosmic strings may be a 
generic prediction of super-string theory [19] and could form in brane collision scenarios 
[20] . In some brane-inflation models like, for instance, the K'^LM^T model [21], D- 
strings are produced as a result of the collision of D-branes. The D-strings are intended 
to be cosmic strings and the D-branes contain a net magnetic flux. The resultant D- 
string network confines that flux in an analogous manner that a vortex network does 
in a superconductor when an external magnetic field is applied. Therefore, a different 
mechanism to KZ must be responsible for the formation of those strings. 

As it was shown in [HI [22] , generally in gauge theories KZ mechanism is only a 
good approximation in those cases where the dynamics of the gauge field is negligible. 
In the opposite limit, to which the examples above correspond, KZ mechanism is not 
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enough and a different explanation is needed for the formation of topological defects. In 
second order phase transitions an explanation in terms of the non-equilibrium dynamics 
of the gauge field has already been found [TH [231 El] • Numerical simulations show the 
formation of clusters of equal sign vortices. 

On the other hand, the possibility for the production of high winding topological 
defects was already suggested in [25]. It was first argued in and then discussed in 
more detail in [22| that they could form as a result of trapping of thermal fluctuations in 
first order phase transitions. We will investigate in this paper the mechanism by which 
thick vortices and heavy cosmic strings indeed form when the underlying symmetry is 
the U{1) gauge symmetry and the phase transition is strong first order. It is in this 
case the joint dynamics of the gauge field together with the complex phase of the order 
parameter - not necessarily out of equilibrium - that gives rise to trapping of magnetic 
flux. 

The reader is warned that throughout this paper we commonly use the term vortex 
to refer to both cosmic strings and vortices regardless of the number of spatial dimensions 
we deal with. Only where the dimensionality is of some relevance we will stress the 
distinction. Also, we will refer to the regions in a minimum energy state either as 
'vacuum regions' or as 'superconducting regions'. Likewise, we will refer to the regions 
in the symmetric phase either as 'symmetric phase regions' or as 'normal phase regions'. 

This paper is organized as follows. Section 2 describes the model and the setup of 
the numerical simulations. In section 3 we present the mechanism by which high-winding 
topological defects form in gauge theories. Section 4 is devoted to the mechanism of 
formation of strips of magnetic flux and chains of single-winding defects. In section 5 
we explain the different properties of the vortex lines in a superconductor film with 
respect to the vortices in the fully two-dimensional theory. In section 6 we argue 
that spontaneous vortices and cosmic strings form in the thermal ensemble. Section 
7 describes some mathematical aspects of the U{1) gauge theory which lead to the 
derivation of Josephson-like effects as well as to the reinterpretation of the so-called 
geodesic rule. 

2. Preliminaries 

2.1. The model 

Let us consider the Ginzburg-Landau (G-L) free energy for superconductors near the 
phase transition point: 



where m is the effective mass of a Cooper pair and /iq is the magnetic permeability of 
free space. 

Let us also consider the Abelian Higgs model, which can be seen as an idealized 




(1) 
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relativistic extension of the above G-L theory. In natural unitsj|| its Lagrangian is 

C = -\f,,F^^ + - Vm- (2) 

Greek indices run from to the number of space dimensions. We will use the indices 
and if: to denote time components and both Latin indices and x, z for the spatial 
components. In the following tensorial formulae only the spatial components apply to 

The charged complex scalar field <^ = |</)| exp {iO) is generically referred to as 'the 
order parameter'. It is interpreted as the quantum wave function of Cooper pairs in 
the G-L theory and as the Higgs field in a particle physics context. The vector field 
is a U{1) gauge field intended either as the electromagnetic field of the usual Maxwell 
theory in the context of superconductivity or as a U(l) field appearing in the spectrum 
of a grand unified theory (not necessarily the photon field). F^j^y is the Maxwell field 
strength and is the covariant derivative operator: 

F^,u = d^A, - d^A,, (3) 
=d^ + ieA^. (4) 

In the G-L theory the value of e is twice the electron charge in agreement with the 
interpretation of a Cooper pair as a two-electron bound state. The Lagrangian ([2]) and 
the free energy ([T]) are invariant under the U{1) gauge transformations 

e ^0 + A, 

A^^A^- (1/e) d^A, (5) 

where A is a real scalar function of space and time. 

In the context of Grand Unified Theories, the model described by the Lagrangian 
([2]) is meant to be just an idealized toy model. When implementing radiative corrections, 
additional couplings of yl^ and (/> to other fields must be considered. Also first order 
time derivatives in the effective equations of motion for both fields, (j) and A^, emerge 
as considering non-zero momenta corrections of order one. 

Concerning dynamical phenomena in superconductivity, the time-dependent 
Ginzburg-Landau equations (TDGL) were derived in [271 [28], [29] cls effective macroscopic 
equations from the microscopic BCS theory. A first order time derivative of emerges 
as a result of the coupling of the electrons in the superconducting phase to the electrons 
in the normal phase. The equation for Ai also contains a first order time derivative as a 
result of the finite conductivity in the normal phase. This is an ohmic current. Thus, the 
dynamics of both fields, and Ai, is approximately diffusive in a superconducting film. 
The TDGL equation for (p does also include a 'small' second order time derivative. The 
equation for Ai must contain a second order time derivative as well. This is necessary 
if the gauge field has to satisfy Maxwell equations out of a superconducting film [21]. 
Therefore, we can explain with analogous equations of motion the dynamics of both a 



I Natural units, c — h — — eo — ks — ^, are used hereafter. 
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superconductor - described by the time-dependent G-L theory - and a particle physics 
phase transition - described by an effective field theory. In the temporal gauge, Aq = 0, 
those equations read 

[du + (Tdt - D,D' + = 0, (6) 

duA., + adtA, + e,,^e\j)^d'A^ = j„ (7) 

where a is the damping rate, eijk is the three dimensional Levi-Civita tensor and 
summation over spatial indices applies if repeated in different location on pairs of tensor 
symbols. The electric current ji = — 2eIm[0*Dj0] is a Noether current which is conserved 
by gauge invariance: 5^[0^(c}^^ + eA^)] = 0. The term —adtAi is an ohmic current jf'^™'. 

Following the original phenomenological Ginzburg-Landau theory in which the free 
energy can be expanded analytically in powers of near the phase transition and the 
fact that many GUTs predict a radiatively corrected effective potential of a similar kind, 
we will work with a scalar potential of the form 

Vm = ml\<P\'-^\M' + M<P\\ (8) 

where mu, k and A are in general positive definite temperature-dependent parameters. 
A first order phase transition takes place at some critical temperature Tc. Below Tc the 
coupling parameters must satisfy the inequalities — 32m^A > 0, — Am'jjX > 
for the existence of a meta-stable minimum at = 0. The true minimum is located 
at \(j)\min = (3k -|- -^/Qk^ — 32Am|^)/8A. In what follows we will refer to |0|mm as the 
vacuum expectation value (v.e.v.) and to the set of absolute minimum energy local 
states as vacuum manifold. Those states are characterized by the conditions 

\4>\ = \4>\min, 

D^cj) = 0. (9) 

The transition takes place by means of spontaneous nucleation of bubbles of true 
vacuum. These bubbles appear as a result of the growth of very special and infrequent 
fluctuations of [201 EI]- Spontaneous symmetry breaking (SSB) occurs as the order 
parameter chooses randomly a point a on the U{1) vacuum manifold such that 
4>vac = \4'\minG^°' (figure [Il^b)). That way, spatial inhomogeneity in the manner the 
symmetry breaks down arises naturally in a first order phase transition. This is a 
requirement for defect formation. 

Nucleation of bubbles and expulsion of magnetic field out of them drive the system 
out of equilibrium. These processes produce latent heat. Restoration of equilibrium 
involves readjustment of the magnetic flux and steady bubble expansion. It is achieved 
by means of the dissipation provided by the damping terms which, at the same time, 
remove the latent heat. 

The damping rates - conductivity in the case of the gauge field - determine the 
time scale of the dynamics of the fields - the steady bubble expansion rate in the case 
of the scalar field. In principle no a priori assumption exists as regarding the damping 
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rates. However they are important phenomenologically since their values determine 
whether the fields evolve in adiabatic, diffusive or dissipative regimes. For the sake of 
simplicity, we will assume that conductivity is uniform in the symmetric (normal) phase 
and we will assign the same damping rate a to both fields. In doing so we guarantee 



2.2. The gauge choice and the gauge-invariant field D^O 

Beyond tree-level, spontaneous breakdown of a local symmetry has been proved to be 
impossible without gauge fixing [32]. In the mean-field theory approximation a similar 
situation holds. Spontaneous symmetry breaking manifests itself in the gauge fixing as 
we will illustrate in section HI 

The fact that topological defects do form implies that a magnetic field gets confined 
around a discrete number of zeros of (j). It results in a multiply-connected vacuum region. 
Along any curve passing through one of the zeros of the complex phase 9 presents a 
discontinuity of value vr. These discontinuities cannot be gauged out because they are 
physical and so gauge-invariant. 

Deep in a simply connected domain of broken phase like a bubble, it is always 
possible to choose the unitary gauge or, more generally, 0{x) = a uniform. In fact, 
because the magnetic fiux is expelled by Meissner effect, it is natural to choose ^4^ = 
there. Conversely, in the symmetric phase region between bubbles where (p = 0, 6 
is ill-defined and so is d^9. As the bubbles of broken phase expand and collide they 
may form a multiply-connected broken phase region. The gauge 6{x) = a can only be 
maintained if this region is not to contain topological defects. If topological defects do 
form, the unitary gauge is not possible because, for any closed loop C around a defect, 
§^ diO dx' ^ holds. 

The above situations illustrate the fact that neither df^6 nor alone are good 
quantities to be tracked in order to explain the formation of topological defects in gauge 
theories for they are not gauge-invariant quantities. In fact, the relevant degrees of 
freedom to be considered in the region where topological defects are to be formed are 
in the gauge-invariant combination of the Goldstone modes and the gauge fielci 



2. 3. Setup of the numerical simulations 

In the numerical simulations presented in this paper we deal mainly with two kinds of 
scenarios. In the first one we show how thick vortices can be obtained when an external 
magnetic fiux gets trapped in a multi-bubble collision. At the same time, chains of 

§ Violation of conservation of charge is not actually a problem in superconductors phenomenology. 
In fact the TDGL equations violate the gauge symmetry explicitly if one considers the system they 
describe as isolated. 

II See section FZ. II for further comments on this definition. 



conservation of charge by fulfilling Gauss' law V ■ E = —q [24j 



a 



(10) 
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single-winding vortices form at the junction of each pair of bubbles. In the second 
scenario we show how chains of vortices and flux strips form in the gap between two 
parallel superconducting blocks in the presence of a net magnetic flux. 

2.3.1. Multi-bubble collision picture 

Nucleation of bubbles is implemented by setting-up disconnected spherical bubbles 
of broken phase. Their initial profile \4>\{r) is approximately that of the bounce solu- 
tion in the thin bubble wall approximation [301 131]. Spontaneous symmetry breaking 
consists of choosing randomly a uniform complex phase value in each bubble: 

= X] — ^ 1 + tanh ^m/,(i?Ar - |r - ri|)j (^cos (a^) + i sin (aj)^ , (11) 
1=1 

where r is the space-coordinate vector, ml is given by fl58l) in section [5|, A*";, is the 
number of bubbles, aj is the initially uniform complex phase of bubble i, stands for 
the position vector of the center of bubble i and is the initial radius of the bubbles. 
Rn is large enough so that bubbles expand. When letting the system evolve with the 
equations of motion, the bubbles grow until colliding. The distance between bubbles 
is large enough so that the bubble-wall profile takes its actual form much before they 
collide. The bubble expansion rate is controlled by the damping rate a appearing in 
(j6]). Soon after bubbles begin to grow the expansion rate v approaches v ~ exp (—a). 

When a uniform magnetic flux is present, nucleation of bubbles requires the 
expulsion of the magnetic flux out of them before they start expanding. We first set up 
a uniform magnetic field B parallel to the y-axis in the gauge = , = ^x. To 
expel flux we let the gauge field evolve according to ([7]) while the scalar field configuration 
is kept fixed at its original form (fTTjl . That way, it is by means of dissipation that the 
system relaxes until getting to the minimum energy configuration compatible with both 
the nucleated bubbles with profile fllip and the presence of a given magnetic flux in 
equilibrium. The optimal damping rate is given by areiax ~ 1/^, where I is the typical 
length of the lattice. The time it takes for the magnetic field to relax is approximately 
treiax ~ ^- After relaxation, both and Ai are evolved according to equations and 
d?]). As it was mentioned above, in order to satisfy conservation of charge, the same 
damping rate a is used in both equations. 

The approach presented here corresponds to the mean-field zero-temperature 
effective theory. However, when a uniform external magnetic field is applied, this field 
resembles a long wave-length magnetic field thermal fluctuation. 

2.3.2. Collision of two- superconducting blocks with magnetic fiux in the junction 



This set-up is meant to mimic the collision between two big bubbles when a mag- 
netic flux gets trapped in the junction. It is a good approximation as long as the radius 
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of the bubbles is much larger than the width of the bubble-wall. The cubic lattice has 
dimensions l^xlyXl^, where < x < l^, < y < ly and < z < l^- The collision axis is 
the X-axis, orthogonal to the walls of the blocks. The gap between the blocks has width 
d and its middle point is localized aX x = 1^/2. A uniform magnetic field parallel to 
the y-axis is confined in the gap. The total flux is quantized with value ^N^- Periodic 
boundary conditions are imposed in the directions z and y and the axial gauge — 
is chosen: 



(p{x,y,z) 



+ 



X 



1 + tanh [mh{lx/2 - d/2 - m^^ 

1 + tanh (mhilx/2 - d/2 - m'^^ 
^cos (q!2 + ^Bd z) + i sin (q!2 -|- eBd z)^ 



— X] 



— \x 



COS (ai) + i sin (ai 



)]( 
-f.l))" 



Ax — Ay — 



{0 for X < lx/2 - d/2 

B{x - lx/2 + d/2) for 1^/2 - d/2 < x < lx/2 + d/2 
Bd for X > lx/2 + d/2, 

where B = '^fjf- and («2 — is the spontaneous gauge-invariant phase difference (see 
below) between both blocks. 

In the fully two dimensional simulations, ly = 0. In the superconductor film model, 
the film thickness /i satisfies /i <^ Ix, ly, Iz- In both models, we refer to the strip between 
blocks of width d along the x-axis, length l^ along the 2;-axis and thickness /i -which 
is zero in the fully two dimensional model- as 'junction'. We will show that it is in 
many aspects similar to a Josephson junction. In the fully three dimensional model, 
as the blocks extend initially all along the y-axis, translational invariance holds at any 
time in the y direction and the problem is actually two dimensional. We refer to the x 
direction as 'transverse direction across the junction' and to the z direction as 'direction 
along the junction'. In the collision of two three-dimensional bubbles along the x axis, 
rotational invariance in the y, z plane rather than |/- invariance holds if no magnetic field 
is applied. We are however interested in the situation where a magnetic field is applied 
in the y direction. Since we take a planar approximation in which the bubbles are big 
enough so that they can be considered as infinite blocks in the z, y plane, the problem 
becomes effectively two dimensional. Therefore, we will maintain the nomenclature we 
use in the fully two dimensional model as referring to the directions along (i.e. the z 
direction) and transverse to (i.e. the x direction) the junction. 

2.3.3. The superconductor film model (SC) 

The superconductor film simulations require extra conditions to account for the finite 
film thickness effects. In the first place, spherical bubbles are replaced with cyhnders of 
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length equal to the thickness of the film /i. More importnatly, vanishing of the variation 
of the free energy ([1]), 6F, with respect to arbitrary variations of the scalar field (p* in 
the superconductor film not only implies ([6]). An additional surface term with the form 

(50*(r)A0(r) ds\r) (12) 



must vanish as well. Q denotes the boundary of the film and ds^ is the differential vector 
normal to the surface at point r. For arbitrary 6(j)*, in order for (fT2l) to vanish, it is 
necessary that the component of Di(j) normal to the film goes to zero at any point of 
the boundary. On the other hand, because the electromagnetic field does not couple to 
any charge particle out of the film, purely Maxwell equations hold in the free space and 
no dissipative term (i.e. a = 0) is to be included there. That is, if the film is placed in 
the a:2;-plane with center at y = and thickness /x, the additional superconductor film 
conditions read: 

=0 for y < —fi/2 or y > fi/2 

a =0 for y < — /u/2 or y > fi/2 

Dy<p = aXy = -/i/2 oi y = /i/2. (13) 

We will be mainly interested in the stability of flux strips and thick vortices. 

The setup for the formation of flux strips is the same as in the previous subsection 
for the confinement of magnetic flux in between two superconducting blocks. 

Formation of thick vortices of winding A^^„ > 1 is achieved by confining initially a 
flux of value "^N^y in a square of side d equal to the typical diameter of the vortex, that 
is d -j^^. If the center of the square is placed at (x = lx/2,y = 0, z = lz/2), the 
initial conditions are: 

(t){x,y,z) = |</)|mm ioi z > IJ2 + d/2 01 z <IJ2- d/2, 



2 



+ 



1 + tanh {rnh{lj,/2 - d/2 - - x)^ 
1 + tanh {mh.{l,j./2 — d/2 — m^^ — |x — 



X (^cos {eBd z) + i sin {eBd z)^ for IJ2 + d/2>z> IJ2 - d/2, 



for 2 > IJ2 + d/2 OT z < IJ2 - d/2 

for X < lx/2 - d/2 and IJ2 + d/2>z> IJ2 -d/2 

A, = I B{x- L/2 + d/2) for IJ2 -d/2<x < k/2 + d/2 

and IJ2 + d/2> z> 1-J2 - d/2 
Ed for X > lx/2 + d/2 and IJ2 + d/2>z> IJ2 - d/2, 
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where B = "^1^2 ■ As the system evolves, a stable, meta-stable or unstable thick vortex 
form. The lattice size is big enough so that finite size effects are negligible. Periodic 
boundary conditions are imposed on the direction y orthogonal to the film. 

Some simulations in this paper differ slightly from the ones presented here. We will 
describe the corresponding setups in each particular case. A description of the lattice 
discretization method is found in Appendix B The values of the parameters used in 
each simulation are compiled in table ICli 



3. Thick vortices and heavy cosmic strings 

3.1. Global Theory 

First suggested by Kibble in [1], the formation of topological defects in a scalar field 
theory with a global continuous symmetry can be explained by means of the spontaneous 
breakdown of the symmetry and the so-called geodesic rule. 

In a first order phase transition of the Abelian global theory, spontaneous symmetry 
breaking takes place as bubbles nucleate at arbitrary points of the U{1) vacuum 
manifold. Since in the global theory any vacuum configuration requires 8^9 = 0, the 
value of the complex phase of the scalar field in each bubble remains uniform as the 
bubble expands. When two bubbles coalesce the above condition cannot be satisfied 
at the contact area unless both bubbles nucleated with the same complex phase value. 
The geodesic rule states that the complex phase 9 of the scalar field must interpolate 
across the junctions between pairs of bubbles in such a way that it traces the shortest 
path in the U{1) manifold. When several bubbles nucleated at different points {ttj} of 
the vacuum manifold meet (figure [I](a)), the complex phase interpolates between each 
(ctj, aj)-pair as shown in figure [U^b). In the case of a three-bubble collision the condition 
for the formation of a single winding topological defect is 

ai + TT < as < 0^2 + vr, ai < 0^2 < «3- (14) 

This condition is fulfilled with probability 1/4. As the complex phase interpolates 
according to the geodesic rule, the gradient energy J dicpd'^cj)* dx across the junctions 
between bubbles, a < x < 6, is minimized: 

pb pb pb 

/ di0<p*dx = 1 di\(t)\d'\<p\ dx + / \(t)\'^di9d'9 dx, (15) 

J a J a J a 

where 9j|0|9*|(;^is the bubble-wall tension and \(j)\'^di9d'^9 is the current self-interaction 
energy densityl^fl More precisely, what is minimized by following the geodesic rule is the 
term associated to di9, that is, the current self-interaction energy density. 

^ The one-dimensional integrations in ()15p are actually either linear or superficial densities of energy 
in two or three spatial dimensions respectively. For brevity, we will refer to them as energy. 



Topological Defects from First Order Gauge Theory Phase Transitions 

V(^ 



11 




Figure 1. (a) Meeting of three bubbles nucleated at different points of the vacuum 
manifold of the global theory, ai,a2,a3. (b) U{1) vacuum manifold. In absence of 
an external magnetic field the complex phase follows the geodesic rule interpolating 
between the pairs ai — a2, as — a2 and ai — in the vacuum manifold. The complex 
phase value 9 is depicted using the colour gradient legend at the top of the figure. 



3.2. Gauge Theory 

In the gauge theory d^6 is not a gauge-invariant quantity anymore. The analogue of 
equation (fT5i) is 



Di(j)D'^* dx 



di\^\d'\(f)\ dx + 



''DSD'd dx. 



(16) 



where the current self-interaction energy density reads e{x) = \(j)\'^Di6D^6. 

When two bubbles collide, a different rule to the geodesic follows: The phase 
gradient D^O and the gauge-invariant phase difference "jab = J^^ D^O dx across the 
junction a < x < b between any pair of bubbles decrease in time until vanishing. In doing 
so, Ai evolves (and so does diO) satisfying the conservation of flux requirement. That 
way, the magnetic flux is expelled, a gradient in the complex phase, dxO, is generated 
and the current self-interaction energy gets to its minimum value. At the same time the 
order parameter takes the v.e.v. across the junction. 

In a multiply-connected broken phase region the above rule applies to any given 
pair of vacuum sectors along each loop C enclosing a simply-connected symmetric phase 
area. Thus, §^ Did dx^ = holds when C is contained in the surrounding vacuum region. 
Continuity in 9 implies diO dx^ = ^nNyj = — §q eAidx^ = —c^e.m.- 

In absence of magnetic field the geodesic rule is a good approximation. In fact it is 
equivalent to the gauge rule above if the gauge Aq = 0, ^^.(t = 0) = is taken, where 
X stands for the direction along the path traversing from one to another bubble. Hence 
the formation of a single vortex between the three bubbles of figure [2] can be explained 
by Kibble mechanism. When the three bubbles of figure [2]^a) coalesce a fiuxon generates 
per bubbles pair in figure [2]^b). If the condition ( |T^ is fulfilled, those fiuxons add up 
constructively at the center of the collision giving rise to a single winding vortex like 
the one in figure EJj^d). Further discussion about the geodesic rule in the gauge theory 
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Figure 2. Fully two-dimensional model simulation. Snapshots of a three-bubble 
collision. Red colour is used for negative values of the magnetic field and green colour 
for positive values. 

in absence of an external magnetic field can be found in section 17.31 

The difference between the gauge rule and the geodesic rule in the gauge theory 
manifests itself in the presence of an external magnetic field. Let us consider for instance 
the phase transition of an infinite superconductor film in an external magnetic field. 
Because the total flux is conserved, the resultant vortex pattern must consist of a lattice 
of vortices of the same sign. On the contrary, the randomness of the complex phase 




Figure 3. (a) Meeting of three bubbles nucleated at different points of the vacuum 
manifold ai,a2,a3 in the presence of an external magnetic flux (in green), (b) Vector 
connection Di0 on the vacuum manifold in the gauge theory. 



pattern obtained by SSB would give rise to the same number of positive and negative 
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vortices according to Kibble's argument. Hereafter we will refer to the above rule in 
gauge theories as geodesic- Meissner rule. 

A typical scenario for the formation of a high winding topological defect is a multi- 
bubble collision in presence of a magnetic flux. As the bubbles of figure H] grow, Di9 
goes to zero inside them and the magnetic flux is expelled. That is the Meissner effect. 
Ultimately, when the bubbles meet and the geodesic-Meissner rule is applied to each 
two-bubble junction, the flux in the symmetric phase region enclosed by the bubbles gets 
trapped and a thick vortex forms at the center of the collision. The main contribution 
to the total flux in the vortex comes from the trapping of part of the magnetic flux 
which is at the moment bubbles nucleate, t^, inside the triangle with contour curve C 
whose vertices are the centers of the bubbles (figure [3](a)). An upper bound is in fact 
given by pa A'^{tN)dxi, where F.T. stands for Flux Trapping. A lower bound is 

^F/T. ~ §c' ^^{'tN)dxi where C is the contour curve of the area enclosed by the bubbles 
when they meet at some time t > tjq (figure [3](a)). The precise amount depends on 
kinematic effects, i.e. on the bubble expansion rate and the dynamics of the magnetic 
field. The upper bound is saturated when either the bubbles expand at a rate 

close to the speed of light or the dynamics of the magnetic field is diffusive. In both 
cases the magnetic flux expelled from the bubbles gets stuck to the bubble walls (within 
a layer of thickness ~ l/o" in the latter case) all the way until they meet. This effect 
appears signaled by the rings in bright green color surrounding the bubbles of figure 
|5l In the opposite limit, when the bubbles expand adiabatically and the magnetic field 
escapes from bubbles at the speed of light, is a better approximation. In general, 
the precise value of $f.t. depends on the bubble expansion rate and the damping rate 
of the magnetic field. In addition, order-one numerical factors are determined by the 
actual geometry of the collision. The quantitative study will be reported in a separate 
article [31]. 

The snapshots of the two-dimensional simulation in figure H] correspond to a soft 
bubble collision in an external magnetic field. After the three bubbles collide gently, i.e. 
f <C 1, a thick vortex gets formed at the very center. Its magnetic flux $e.m. lies in the 
range < ^^e.m. < because a is chosen so that the dynamics of |(^| is close to 
adiabatic while the dynamics of the magnetic field is diffusive. It is easy to calculate 
its value by simply counting the times the complex phase increases 27r units around the 
vortex. 

The contribution of Kibble mechanism through SSB is negligible in the presence of 
an intense magnetic field, that is, for ^f.t. ^ 27r/e. Let us consider a fixed value for 
the initial magnetic field in the setup of figure |4](a). If one chooses randomly the initial 
{9i} values of the complex phases of the three bubbles the winding number of the vortex 
that gets formed at the center of the collision takes just one of the three possible values 
■^^F.T. ± 1- Only if the contributions due to Kibble mechanism and Flux Trapping were 
additive, according to the probability for condition f|T^ to be satisfied, the probability 
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Figure 4. Fully two-dimensional simulation. Snapshots of a three-bubble collision 
in the presence of an external magnetic flux (in green). The colours in the bubbles 
stand for the complex phase value of (j> according to the colour legend. In (d), a thick 
winding-five vortex gets formed at the center of the collision. 



P of getting each value would be: 

Pi^'^F.T.) = 3/4, 
Pi^^F.T. + 1) = 1/8, 

Pi^^F.T. - 1) = 1/8. (17) 

However we will show at the end of section H] that Kibble mechanism and Flux Trapping 
are not actually additive. 

If the coUision is violent, i.e. f ~ 1, as it is the case in figure 0, the flux trapped 
corresponds roughly to The non-equilibrium dynamics is relevant, reason why 

the flux trapped is not only confined in the form of a thick vortex at the very center of 
the collision. On the contrary, chains of single-winding vortices form at the junctions. 
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Figure 5. Fully three-dimensional model simulation. Snapshots of the cross-section 
of a violent three-bubble collision in the presence of an external magnetic flux. The 
cross-section is taken on the collision plane. Information about the complex phase is 
omitted. The broken phase bubbles appear in black while the magnetic flux density 
is depicted in green colour gradient. A thick cosmic string form at the center while 
chains of single winding strings form at the junctions. 

4. Magnetic flux strips and single vortex chains 

By causality, when two bubbles collide in an external magnetic field, a finite amount of 
magnetic fiux gets trapped at the junction. Let us assume in first approximation that 
small structure effects -i.e. effects at scales comparable to the bubble wall thickness 
~ m'^^- can be neglected in the dynamics of the bubbles. That way the bubbles expand 
isotropically at constant rate v at all times, even after their walls overlap. Let us assume 
the less favorable situation for magnetic flux to get trapped. That is, let us assume that 
the magnetic flux consists of massless magnetic particles which are expelled from the 
surface of the bubbles at the speed of light as if they were light-rays. That way we are 
neglecting several effects. First, we are assuming no interaction between the magnetic 
field and the scalar field at small scales of the order of ~ f^h^- Second, we are ignoring 
the long range magnetic field interaction, which we will show later on is present in the 
superconducting model. Third, we are neglecting any diffraction effect. All these effects 
contribute positively to the trapping of magnetic flux otherwise. As a consequence, the 
estimate of the flux trapped that we proceed to calculate disregarding these effects is a 
lower bound. 
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Figure 6. Trapping of a 'magnetic ray' by causality. A 'magnetic particle' is expelled 
at the speed of light c from the surface of the left bubble at {xi{to), zi{to)). It gets 
trapped by causality at some point (x = D/2,z) along the junction between the two 
bubbles. The centers of the bubbles are separated a distance D and both bubbles 
expand at constant rate v. 

and expanding at constant rate v as shown in figure El A magnetic particle abandons 
the surface of the bubble on the left at height -21(^0) from the collision axis when the 
bubble radius is tq = vt^. The ray experiences several reflections as it strikes several 
times on each bubble. These reflections may, or may not, follow Snell's law. In the case 
they do, the process is referred to in [M] as magnetic optic approach. Nevertheless, in 
the worst of the cases the exit angle is arbitrarily large at each reflection. However, 
the distance the magnetic particle travels in a time interval At = t — t^ > cannot 
exceed c(t — to)- If t > -D/2f, bubbles meet and overlap during that time interval such 
that the vertex of the overlap is located at height Z2^'^*'^^ {t) . Because the speed of the 
vertex along the vertical axis z is initially superluminal, it is possible that for some time 
interval T/-, 

c(t/ -to) < ||(xi(to),^i(to)) - P/2,^r*'^(t^))||, tf e Tf, (18) 

where (xi(to), ^i(to)) stands for the coordinates of the point from which the ray is first 
expelled and (D/2, 2;2^''*^^(t/)) stands for the coordinates of the vertex. If inequality 
ffTSjl is satisfied for some tf E T then the light-ray must have got trapped at some point 
of the junction {x = 0) where bubbles overlap. 
The vertex height, Zg*^*"**^^, evolves in time as 



^vertex (^^^ = ^vH^-D^i, t > D/2v. 



(19) 
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The maximum height Z2 at which a photon expelled at (xi(to), -2i(^o)) can get along the 
junction at time t > to is 

Z2{t)\i,,ito),z,(to)) = zi{to) + v/c2(t-to)2-(xi(to)-/^/2)2, t>to + (l/c)(D/2 - Xi(to)). (20) 

From fl20|) one reads that the rate at which Z2 grows is initially infinite. However, from 
the analysis of its second derivative, it tends asymptotically to c faster than the first 
derivative of (fT9|) tends to v. As c > v, inequality ( |T8l) may get saturated at time t^""^ 
for some pair {xi{to), zf'^"-^(tQ)) such that 



Z2{tf "'')|(xi(to),^f--(to)) 



Therefore, the fiux in the differential area zf^"'^{xi)dx gets trapped. We have solved the 
equations above numerically discretizing the x-axis with space step 6x = D/2N. We 
search for zf'^"'^ for each x\ = ix 5x, i G [0, A^]. Starting with a uniform initial magnetic 
field Bq, the fiux trapped by causality is given by 

N 

^Caus.trapped ^ ^ zf^^^ix^Sx. (22) 

i=0 

Figure [7] shows the causality bound estimate (curve (b)) in comparison to the actual data 
(curve (a)) of the simulation of a two bubble collision in an initially uniform magnetic 
field Bq. Quantization of the magnetic fiux has been taken into account in the causality 
bound computation. A damping rate a is used to slow down the bubble expansion rate 
V through ([6]). However, no damping term is included in the equation for Aj. That way 
the dynamics of Ai is neither dissipative nor diffusive. 

Having shown that magnetic flux gets necessarily trapped at the junction by 
kinematical reasons (modulo quantization effects), we aim to explain how it gets confined 
in the form of vortex chains like those in figure [5](d). This involves the study of small 
structure effects, that is, the dynamics of |0| and Di9 at the junctions. 

Let us simplify the problem by considering that the collision area of two big bubbles 
can be approximated by two parallel superconducting blocks separated by a gap of width 
d in the transverse direction x in which a net magnetic flux is placed. For the sake 
of simplicity we will work with the fully two-dimensional theory. Periodic boundary 
conditions (pbc) are imposed in the direction z along the junction. The total flux in the 
gap is quantized. Its value is ^N^, where is an integer. The reason for imposing pbc 
is to conflne the flux as it would happen when two big bubbles coalesce (by dynamical 
reasons in the latter case). 

Given two disconnected broken phase blocks we introduce a uniform magnetic flux 
in between. With no loss of generality we will choose the gauge A^, = at t = as 
shown in flgureO In this gauge the fleld Ai is continuous and so are di6 and Di6 in the 



vertex fj.Alax\ 
^2 V^f ) ) 

J ^vertex /4.\ 

dt ^21) 




Figure 7. Number of single-winding vortices trapped in a two-bubble collision. The 
bubbles expand at constant rate v. Initially, there exists a uniform magnetic field 
B = 7.76 • 10^^/e and the distance between the centers of the bubbles is D = 300 
in lattice spacing units. Curve (a) corresponds to the actual fully two-dimensional 
simulation. Curve (b) corresponds to the causality bound estimate. 




Figure 8. Two superconducting blocks separated by a gap of symmetric phase. The 
magnetic flux in the gap (in green) is ^N^, with iV^ = 3. The gauge is chosen in 
such a way that Di — in the superconducting blocks. The complex phase value in 
the gauge below is depicted according to the color gradient legend. 
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areas of broken phase. Furthermore, the gauge-invariant phase difference 'jabi^) across 
the junction when the blocks meet at time t = is equivalent to the complex phase 
difference A6ab{z). That is 

1ab{z) = / DJ{z) dx 

J a 

= / dx9{z) dx = A9ab{z) in the gauge = 0. (23) 

J a 

Therefore, according to the gauge choice in figure [HI 

lab{z) = {as - aA) - (N^/R) z, (24) 

where (a^ — a^) is the spontaneous phase difference and —{N^/R) z is the induced 
phase difference)^ After fixing the gauge there exists still the freedom in choosing 
Aoafe = {as — OiA)- This freedom reflects the invariance of the system along the 
junction (z-invariance) when the blocks are separated. The equality 7afe = Actab 
holds in absence of magnetic field. Thus, the arbitrary choice of Auab corresponds 
to the spontaneous breakdown of the U{1) symmetry as announced in section [2l2l The 
spontaneous difference Aaab is physical and manifests itself as blocks overlap and the 
2;-invariance gets broken. At the moment the blocks get in contact, the physical local 
gauge-invariant quantities are not translation invariant anymore along the junction due 



to the presence of the magnetic field (see Appendix A for a proof) 



The scalar field profile falls off exponentially within a typical length equal to the 
inverse of its mass l/ruh. Consequently, it is necessary that the size of the gap between 
the blocks d is less than dg ~ 2/mh so that the blocks overlap and interact with each 
other. 

If d ^ dg, the blocks behave as fronts of big bubbles whose expansion is driven by 
the potential V^(|0|) and the bubble- wall tension The 2;-invariance of |0| and 

B along the junction persists and no current traverses the gap. Consequently no vortex 
can form. On the other hand the trapped flux gives rise to a repulsive force which, if 
d > dg, can balance the expansion of the blocks towards each other. If one neglects the 
bubble-wall tension in first approximation it is easy to calculate the critical value for 
the magnetic field necessary to prevent the blocks from approaching each other beyond 
some distance greater than dg. At zero-temperature the equilibrium configuration 
corresponds to the minimum energy configuration compatible with a net magnetic flux 
in the junction. In the fully two-dimensional theory the magnetic field B does not 
present long-range interaction in the symmetric phase and the equilibrium configuration 
corresponds to a magnetic flux density uniformly distributed in the junction. Therefore, 
for a fixed value of the magnetic flux, the magnetic energy depends only on the size of 
the junction. The equilibrium distance deq between the two superconducting blocks of 
figure [8] is the distance at which the expansion of the blocks is compensated by the 

See section FZ. II for further explanation about this decomposition. 
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magnetic pressure, 

, _ ^E.M. 

where 27tR is the length of the blocks along the junction, ^e.m. is the total flux and 
Vmin stands for |V(|0|mm)|- Requiring deq = dg we calculate the minimum value for the 
magnetic flux so that the blocks never overlap and the flux strip remains straight, that 
is, 

47r 



^T^''^''' = —R^W~. (26) 

If the value of $ em. in the gap is less than <|)^^^'™" the small scale structure effects 
are relevant and the dynamics of both |0| and DiQ cooperate towards the formation of 
topological defects. A detailed analysis of the generation of Josephson-like currents can 
be found in [33]. In what follows we just make use of the essential results there. 

At the moment the blocks get in contact, ^ah{z) = /^^/2 ^^^i^i z) dx tt for some 
z, where x is the typical length of the region where D^O ^ and dx\<f)\ ~ 0, x <^ (i. 
This implies that Dx6{z) is a function of "yabiz) and the equations of motion incorporate 
non-local effects and non-linear terms in D^O. In particular, a Josephson-like coupling 
shows up (see also sections 17.21 and 17731) . The quantity to be studied is 'Jabiz) which, 
right at the moment the blocks get in contact, follows a perturbed sine-Gordon equation 
in 1 + 1 dimensions 



dttlab - dzzlab + crdaab + 2e^|0|o sin7af, = 0, (27) 

where |(/)|o is the independent factor of the norm of the scalar field at the time the 
blocks meet, i.e. |0|^2,|<^/2 — I'i^lo ^ fi^)- I0lo i^ referred to as 'reduced Cooper 

pairs density'. For x <C rf, it is 

|</)|o~(l/2)|0|™„e-'^/l (28) 

Therefore, initially, the vortex formation dynamics corresponds to that of one- 
dimensional sine-Gordon solitons. At this stage the current self-interaction energy 
density depends on z as 

^|0|^(l-cosM;.)]) =1 

X X 
At \(f)\o fixed, equation implies that Sab gets maximum at those points Zn satisfying 
1ab{zn) = i^2n + 1)tt, which are 



eM) ~ -\<P\l{^-cos[-f,t{z)]) = -\<l)\lsm'[^ab{z)/2]. (29) 



R 

Zn = —[ae - aA - i2n + 1)tt], < n < N^, (30) 

where the linear flux density is expressed in terms of the winding number, '^'^^'J^' = 

When the blocks get in contact, 4>\x<o and (f)\x>o interfere at the middle of the gap 
X = giving rise to an interference pattern along the junction as 

|0|Lo ^ 2|0|2(l + cos [jabiz)]) = mlcos' [jab{z)/2]. (31) 
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Such an interference is constructive at the points Zn' such that ■yabi^n') = 27rn' (equal 
colours according to the color gradient legend in figure [8]). On the contrary it is totally 
destructive if ( l30l) is satisfied. As a result A^^ minima of |0| turn up at the points 
{{x = 0,z = Zn)} when the blocks first meet, (see figure lAlT b)). 

In [33] it is shown that \(f)\o is not z- invariant anymore either as |0| follows the 
equation of motion fl88l) . We will just use here a minimum energy principle to argue for 
the formation of vortices. In view of formula fl3T|) . the potential fl29|) gets minimum at 
the same points where |0| is maximum and vice verse. The points {{x = 0,z = Zn)} are 
in fact centers of the orbits of the current fiow (figure [ATT d)) and attractors of magnetic 
fiux. The current flow ~ DO turns around the minima of |0| while the magnetic field, 
B = (1/e) V X DO, goes straight into them. They are seeds for the formation of a chain of 
single winding vortices if an additional condition holds (see below) . Equation ( !27l) must 
be intended as a good starting point to evaluate the initial evolution of ■jab- Afterwards, 
additional non-linear terms enter and the problem becomes two-dimensional in space. 

It might still happen that the magnetic repulsion stops the vortex formation process. 
Unit-winding vortices are rotationally symmetric with a radius that scales with the 
London penetration length which is equivalent to the inverse of the photon mass |36j : 
^ l/nij = (-\/2e|0|rrtm)~^. Therefore the inequality tcR > A'^r„ must be satisfied in 
order to accommodate Ny^ single vortices in a chain of length 27tR. This condition 
imposes an additional constrain on the maximum value of the fiux trapped for a chain 
of single-winding vortices to form: 

^chain,max ^TT R r,3/2 2T3\i\ fQO\ 
^E.M. = ^ 2 ' TT R\(p\min- (32) 

6 

In the intermediate situation the magnetic fiux lies in the range > > 

^cham,max ^ In that casc the minima oi\(j)\ do form but single vortices do not. Instead, 
periodic fiux 'ripples' appear along the junction between blocks. Figure [9] shows the 
three typical cases. The stability/meta-stability of the configurations [9](b) and Ws) 
depends on the dimensionality of the theory. That is, fiux 'ripples' and chains of single- 
winding vortices are not stable against condensation into thicker vortices or heavier 
cosmic strings. However they can be meta-stable for high enough fiux density. That 
meta-stability is stronger in the superconducting model where a long-range magnetic 
interaction is present (see next section). Figure [TUT a) shows a membrane of magnetic 
fiux in between two superconducting blocks. The fiux membrane contains nine winding 
units of fiux. They manifest in the form of 'ripples' in figure [TOT b). The equilibrium 
of such configuration is just maintained by the periodicity in the direction z along the 
junction. However, when we introduce a small perturbation <^ \4>\min at z = 0, the 
periodicity is broken and the structure becomes unstable. The vortices gather in groups 
of winding three (figure [TOlfc)) and three thick vortices of winding three get formed 
(figure [lOKd)). 

The spontaneous phase difference as — a a is indeed physical. It gives rise to 
a shift in the location of the vortices along the junction according to ^2M . In a two- 
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(a) X (b) (c) 



Figure 9. Fully two-dimensional model simulation. Typical magnetic flux structures 
at the junction between two superconducting blocks. The lengths of the junctions 
27ri?() increase from left to right while the magnetic flux is 67r/e and the spontaneous 
phase difference is 37r/4 in all the cases, (a) Flux strip, (b) Flux 'ripples', (c) Vortex 
chain. 




Figure 10. Superconducting model simulation. Snapshots of the evolution of a 
membrane {light grey) of magnetic flux 187r/e between two blocks of broken phase 
[dark grey). Periodic boundary conditions are imposed in the directions y, z. The film 
thickness is /i = 3 in lattice spacing units. Light grey is used to depict isosurfaces of 
equal magnetic energy density. Breakdown of periodicity along the z-axis between (b) 
and (c) leads to the formation of three vortices of winding three in (d). 
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bubble collision, whether a vortex forms at either side with respect to the collision axis 
depends on this shift. In figure [TT] two bubbles 'sandwich' a strip of magnetic flux as 
they coalesce. The relative shift difference of value Sir /A between the cases (a) — (6) and 
(c) — (d) gives rise to the production of one more vortex along the junction in figure 
ITTlfb) with respect to figure fTTTd). This is the main effect that SSB has in the total 




Figure 11. Fully two-dimensional model simulation, (a) A strip of uniform magnetic 
flux 367r/e is placed between two bubbles with spontaneous phase difference Aai — 0. 
(b) After the collision of the bubbles in (a) a chain of five single- winding vortices forms 
at the junction, (c) A strip of uniform magnetic flux 367r/e is placed between two 
bubbles with spontaneous phase difference Aa2 — 3tt/A. (d) After the collision of the 
bubbles in (c) a chain of four single-winding vortices forms at the junction. 

flux trapped in the thick vortex which forms in a multi-bubble collision. The fact that 
a single-winding piece ends up either forming part of the thick vortex at the center or 
being expelled depends on this shift. 

There is also another way in which SSB and flux trapping mechanism interact. 
Let us consider the current traversing the gap between the two bubbles in the presence 
of a magnetic flux. The total current follows an interference pattern (see formula fl99|) . 
section [7]2]) and oscillates as a function of the total flux. Its amplitude decreases notably 
as the flux increases and it is zero if the total flux in the junction takes integral values of 
27r/e. As a result of that current a fluxon ring would develop around the collision axis 
giving rise to an extra contribution to the magnetic flux in the thick vortex (similar to 
what happens in figure [2] in absence of magnetic field). However, because that current is 
strongly suppressed for high values of the magnetic flux in the junction, we expect that 
fluxon to be very weak. Consequently its contribution to the flux in the thick vortex is 
negligible. For instance, in the collision of the three bubbles of figure HI the probability 
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for the final winding in the thick vortex to be ^^f.t. must be greater than the value 
3/4 in f|T71) due to this interference effect. We conclude that he contributions from flux 
trapping and KZ mechanism are not additive. 



5. The superconductor film model 

The model which describes a superconducting film differs from the fully two-dimensional 
G-L model in the fact that the order parameter (p is confined in a thin film while the 
electromagnetic potential Ai lives in the three dimensional bulk space. The magnetic 
field is not a pseudoscalar but a three-component pseudovector. It propagates in the 
bulk and induces a long-range interaction on the magnetic flux density in the film which 
is absent in the fully two-dimensional case [23]. 

Let us take a thin superconductor film of thickness /i, placed aX y = and 
parallel to the xz-plane. If a net magnetic flux is present, it gets confined in areas 
of symmetric phase. For a stationary configuration, the electromagnetic energy reduces 
to the magnetic energy confined in a volume V whose cross section at ?/ = is the 
symmetric phase area on the film, E^^-J^^- = | /y d^r BiB\ As E^^Jl^- is quadratic in Bi, 
its most general form in terms of the film magnetic flux density By{r) reads 

Efd'- = lf I d'rd'f^ By{r)V^^,.{r,f^)By{f^\ (33) 

J sym. J sym. 

where two-dimensional vectors with components x and z in the film are labeled with 
arrows on the top, Vmag.if.r') is the induced magnetostatic potential interaction and 
the subindex sym. stands for 'symmetric phase area'. The stationary configuration 
corresponds to that for which fl33|) is minimum. That is, By{f) such that 



6By{r) 







d^f'Vmag.{ry)By{n, (34) 

' sym. 

subject to the presence of a fixed magnetic flux ^e.m. = /^ym d'^r By{f). We first 
calculate Vmag.i^,!^) from the stationary Ampere's law. 



(V X B)i = ji. (35) 

Because the current lives in the superconducting film on the plane y = 0, only the 
components x,z are non-zero in the vector equation fl5Sl) . On the other hand, because 
the magnetic field is orthogonal to the film at y = and the flux gets trapped in the 
symmetric phase regions. Ampere's law reads 

[VxB{r)], = [VxBy{f)]„ reUy=l, z = x,z, (36) 

where stands for 'symmetric phase area in the plane y = 0\ 

As a first approximation we will consider the film infinitely thin, that is, with zero 
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thickness, fi = 0. That way, the ?/-component of the magnetic field on the film can be 
written in Fourier space as 

B^=%k) = j £r By{r)5{ry)exp{ik-r), = 0. (37) 

It is more convenient to write fl36|) in function of the gauge field Aj(r), 

[V X V X A{r)], = [V X By{r}],, f G n^=l, t = x,z. (38) 
In the Coulomb gauge, VjA* = 0, 

- VM,(r) = [V X 5,(r)]„ feuy=l, i = x,z. (39) 
We solve for Ai{r) using the Green's function of the operator V^, that is, Gj^ such that 

VlGl\r,r')=5,,5^^\r-r'). (40) 

In Fourier space, 

Gl\k) = %. (41) 

The gauge field can be computed as 

Mr) = - f dr'y6{r'y) f rfV-'Gj'(r, r')[V^' x Byif')]^ (42) 

J J sym. 

and likewise the magnetic field, making use of ( 14T1) and expanding By{r) in Fourier 
modes, 

B,(r) = [V X A(r% 



where Sijk stands for the Levi-Civita tensor, the subscript y is fixed, m, n take labels 
{x,z} and i,j,k run over {x,y,z}. Next, we define the 'pseudo-Green's function' 

g.(r, n^-j dr;S(r') J ^||e"-<'-")e'--.,„4»„^ (44) 



and its corresponding Fourier transform 



^f(fc,p) = -e,,^e-^, (45) 



such that 

Bi{r) = [ d^7'g^{ry)By{7'). (46) 
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E 



E.M. 

stat. 



- I d^rBiB' 

1 
2 



c/V / c/V / dh By{f^)gf {r,r')Gy\r,r")By{f^'), (47) 

I. J sym. J 

where f denotes adjoint. From fl33|l . one reads immediately 



Vr 



•mag.y i 



r ,r ] 



d'rgf{r,r)gy\r,r' 



(4^ 



Again, we make use of the Fourier transforms to write, in the infinitely thin film limit 
/i ^ 0, 



d^p d'^s 



By{p) / By{S 



d^k d^q 
(27r)3 (27r)3 



X gfik,p)gy\q,s) 



(27r)2 (27r)2 
/i ^ 0. 



(49) 



After the integration of s (which yields a product of delta functions) we get to 



E. 



E.M. 
stat. 



- / d^fi"d^f2' 



d^k ~ 



By{k)By{-k)e''- 



ik-{fi"-ri") 



(50) 



(27r)3 -Uf-^aft-ym-j/d ^4 

where k'^ is the 4*'' power of the norm of the three-component wave vector k; r, b, m, d 
take labels {x,z} and i,j,a run over {x,y,z}. By applying the identity etjke''''"^ = 
^i^T ~ ^T^j integrating the mode ky , it is easy to get to 
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Comparing (ISTj) with (14711481) . we identify 

d^k 



(2vr)2^ 4'^'^ 
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(51) 

(52) 
(53) 



It is remarkable that VJi^^g {7\ ) equals the inverse of the two-point function 



{By[fi')By{f*2)) in the high temperature field theory, normalized to ksT = 1 [23]. The 
reason being that the long wavelength spectrum of fluctuations of the gauge field at 
high temperature -that is, the Ray leigh- Jeans spectrum of balckbody radiation- equals 
the Green's function of the stationary Ampere's equation ( l39l) in the symmetric phase. 

Should we take into account the finiteness of the thickness of the film, /i, on the 
magnetic flux density By{k), we would integrate those modes for which ky < 1/yU. We 
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may use a Gaussian function to suppress shorter wavelength modes. Replacing the delta 
function S{ry) in the integrand of fl37|) and thereafter with the Gaussian distribution 
exp (— r^/2yU^), we define the film magnetic flux density as 



2lTfl 



B''>^{k) = [ d^r exp {-rl/2fi^)By{r) exp {ik ■ f) 

J V2vru 



'2TTfI 

This way, analog of ( l52|) with fi > reads 



(54) 



V^>J.irr, rT) = I ^e^^^-(-"--")[iErfc(/.|^|) \k\]-\ (55) 



where Erfc stands for 'complementary error function'. If the typical length of the 
symmetric phase area is much larger than the film thickness, the long wavelength limit 
\k\ <C 1/fi can be taken in good approximation such that Erfc(//|/c|) ~ 1 and the 
approximate expression (152|1 is recovered. 

We end up with the result that, in a superconductor film, vortices present a 
long-range interaction associated to their magnetic flux. Also, the Coulombian-like 
interaction K«ag.(^;^) implies that the magnetic flux density on the symmetric phase 
areas of the film behaves as if it were an electric charge density on the surface of an 
ordinary metal conductor. 

When a magnetic flux is trapped in a symmetric phase region, the stationary 
configuration corresponds to that for which ( |33l) is minimum subject to the existence of 
a total magnetic flux ^e.al- That is, By{r) such that 



5By{r) J J sym. 




d^rcPr 



\B,(f)V^,.(r. r-)B,(f) H- H S^f- r-)B,(f) = 0. (56) 



where if is a lagrangian multiplier and Vmag.{r, r') is given by (!53l) . The solution of (!56ll 
is 

By (f) = ~H I rfV-' (f , r^) , (57) 

J sym. 

with the constraint ^e.m. = J J^y^ d'^fd'^r' 6^{f— r')By{r). 
5.1. Stability and meta- stability of thick vortices 

For the purpose of studying the dynamics of vortices, there are in principle two length 
scales to be considered, the coherence length ^ and the London penetration depth A^. 
^ is equal to the inverse of the mass of the order parameter at |0|mj„, that is, ^ = l/nih, 
where 

™N4fl^»)|uH*U.. (58) 



For the potential in ([8]), m\ = mj^ — 3fi;|0|mm + 6A|0|^j„. Notice that even though 
ruh is of the order of mn they differ by an order one numerical factor. is the 
inverse of the photon mass, Xl = = (A/2e|0|mm)^^- The Ginzburg parameter. 
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= ^h/i = 'mhl'i^-yi depends on both the value of |0|mm and the curvature of the 
effective potential rrih at |<^|mm- Finite temperature radiative corrections [37] and the 
presence of an external magnetic field [38] may give rise to variations in kq- 

In a type-I superconductor, where < ^, there exists an attractive vortex-vortex 
force at distances shorter than ^ associated to the order parameter. At distances shorter 
than \i there is also a magnetic force interaction which is repulsive for equal sign 
vortices. If no other force exists, two vortices of the same sign separated by a distance 
less than l/m^ can gather to give rise to a stable thick vortex. That is the case in the 
fully two and three-dimensional theories [39] . 

In a real type-I superconductor film, the induced long-range repulsive force (l53ll 
between either equal sign vortices or equal sign flux densities, modifies the picture above. 
In a thick film, 3> ^, the long-range Coulombian interaction is dominant at very long 
distances. Therefore, a meta-stable cluster of vortices of the same sign form. If the film 
thickness is of the order of ^ and still greater than A^,, it happens that high winding 
vortices are stable up to some critical value of the winding number. For higher winding 
they can be meta-stable . Finally, if the film is thin, fi ^ Xl, the Coulombian repulsive 
interaction acts at distances r > Xl and so only single- winding vortices are stable and 
an Abrikosov lattice form in an external magnetic field . The last case corresponds to 
the situation in which a thin type-I bulk superconductor film behaves as type-II in a 
perpendicular external magnetic field [^ 1^ fi2] . 

Let us concentrate on the case where both stable and meta-stable thick vortices can 
form. The winding number of the stable vortices cannot be arbitrarily high. Beyond 
some critical winding value, they break up into smaller pieces. The argument for this 
to happen is as follows. The total energy of a stationary winding N^j vortex contains 
three terms, 

E{N^, r,) = C{fi/r,)N^ r;' + fiV rl + T r„ (59) 

where C^^/r^) is a geometrical factor which depends on the exact form of Kraag.(r, 
through /i/r„ - formula (|53|) is only exact in the thin film limit /i/r„ ^ -, = nVmin 
and T stands for 2it times the tension of the vortex wall. Stability requires 

dE dC 

— \,sta,. = -—^Nlrf - C Nlrf + 2fiVr, + T = 0. (60) 

Considering the film thickness fixed and the vortex radius large in comparison to both 
fj, and Xl so that formula (153|) is a good approximation and the energy associated to the 
vortex wall tension is negligible (i.e. Tr„ -C CN^ r~^, fiVr^), the first and fourth terms 
in (|60|1 can be ignored and 
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stab. ^ ^ Kf'^/'^ 
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E'^^'iN^) ^ (2^/3 ^ 2-2/3)^2/3^1/3^1/3 (61) 

Consequently, 

E'^^^iN^) + E'^^^il) < E'^^^iN^ + 1) V A^^ (62) 
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and multiple winding vortices are not stable. 

Let us illustrate this phenomenon by analyzing the stationary configuration of the 
magnetic flux density in a thick vortex. For sufficiently large values of A^^, the vortex 
radius scales as r„ oc A^^^^ and the magnetic energy is dominated by the long-range 
interaction of the flux density in (|53ll . In that case the minimum-energy equilibrium 
conflguration of the magnetic flux ^e.m. = in a vortex approaches the distribution 

of an electric charge placed on an ordinary conducting surface of circular shape and 
radius [13]. Solving for By{r) in flFTl) . 



where r is the distance to the center of the disc. Formula (|63l) gets maximum along 
the periphery of the circle rather than at the center as would be the case in a fully- 
two dimensional vortex. As A^^ is increased, Vy goes to inflnity while By\r=o tends to 
zero. Such a tendency maximizes the interaction energy between the magnetic fleld and 
the scalar particles through (flGl) . The magnetic flux density By goes to zero at r = 
where is zero while By increases as r goes to where |0|r=r„ ~ |0| mm, making the 
configuration energetically unfavourable. 

However, given a pair of values C{fi/ry), /i in fl60l) . the energy associated to the 
vortex wall tension may compensate the energy difference in the inequality ( l62i) for small 
values of A^^ and r^. If that is the case, inequality (162!) can be inverted. Figures [12] and 
[T3] show this phenomenon in numerical simulations. In figure [12] the magnetic fiux is 
Gn/e and A^, is approximately one in lattice spacing units. When the superconducting 
film is thick (/x = 21 in lattice spacing units) the magnetic fiux profile presents a 
maximum at the very center and falls off exponentially towards the periphery of the 
vortex. As the thickness decreases, the profile of |0| does not change sensibly. However 
the magnetic fiux profile does approach the form of formula (1631) and the maximum 
magnetic fiux density migrates towards the periphery of the vortex. Finally, for a 
thickness ^ = 1 ^ Xl the winding-three vortex splits up in three single-winding pieces. 
In figure [12] only the profile of one of them is plotted. While the winding-three vortex 
is stable for fj, = 21, it is meta-stable for /i = 9, 5, 3 and unstable for /i = 1. Figure 
[13] shows how the energy per winding number of thick vortices E{Nw)/Nu, scales with 
the winding number A^^„ in comparison to the energy of a single- winding vortex E{1). 
For /i = 1 only single- windings are stable. For 15 > ;U > 1 only single- windings are 
still stable but numerical simulations show that winding two and three vortices are 
meta-stable. It is difficult to test how much meta-stable higher windings are because 
the number of possible combinations and geometries in which a thick vortex can decay 
increases quickly with its winding. For yU = 15 a winding-two vortex is stable against 
decay into single-winding pieces. For /i = 21, winding- four vortices are stable and, at 
least up to A^^ = 6, thick vortices do not split up in single- winding pieces. A winding-six 
is however meta-stable and can decay into a pair of stable winding-three vortices. 

We expect that, by the mechanism explained in previous sections, a variety 
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for r > r, 
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Figure 12. Superconducting film simulation, (a) Profile of the scalar field, (b) Profile 

of the magnetic flux density at the meridional plane of the film y = 0. ^e.m. = 67r/e 
for film thickness /x = 3, 5, 9, 21. $e.m. = 2n/e for the curve = 1. 
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Figure 13. Superconducting film simulation. Scaling of the energy per winding 
number of thick vortices with the winding number. E(\) stands for the energy of a 
single-winding vortex, [i is the film thickness in lattice spacing units ^x. Because in 
our simulations Ai, « /x- values are approximately in units of Ai,. The dashed line 
marks stability against decay into single- winding pieces. Vortex stability is guaranteed 
at points of non-positive slope. 
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of stable and meta-stable vortices form during the phase transition of a type-I 
superconductor film in an external magnetic field. Because of the inhomogeneous nature 
of a first order phase transition, the primordial lattice can be very different to the one 
corresponding to the minimum energy configuration. Its evolution in the first stages 
must be dominated by the long-range magnetic interaction. It is also very likely that, 
by the enormous amount of meta-stable intermediate states, the evolution of the vortex 
lattice gets stuck to some meta-stable configuration. It is only for very thin films that the 
only stable vortices are single-winding and the final configuration is a stable Abrikosov 
lattice. The Abrikosov lattice configuration is unique up to the fixing of its center of 
mass. The center of mass consists of two degrees of freedom which correspond to the 
two-dimensional analogue of the spontaneous phase difference Aaab that determines 
the position of the single-winding vortices in a finite vortex chain. Thus, they are also 
vestiges of the spontaneous breakdown of the gauge symmetry. 



5.2. M eta- stability of flux strips 

The strips of magnetic flux described in section H] are generally unstable in the fully two 
and three-dimensional models. For thick superconducting films like the one in figure 
[To] where /i ~ BA^ they may be unstable as well. They can however be meta-stable 
for sufficiently thin superconducting films and strong enough magnetic field. That is so 
even if their width is much less than 2/m/i and short distance effects become relevant. 
Numerical simulations show that, for a fixed value of the linear magnetic flux density 
higher than $^'^"'™"^/27ri? (2tiR is the length of the strip), the resultant strip is meta- 
stable for a film thickness less than some critical value. Conversely, given a fixed value for 
the film thickness, there exists a critical magnetic flux linear density for which the strip is 
meta-stable. As an example let us take the same setup of figure [TO] for a superconductor 
film of finite thickness /i. Using the parameter values in Appendix C if /i is set to 3 in 
lattice spacing units, the resultant strip is unstable for = 9 and three thick vortices 
of winding three form as in figure [THJ However, if we increase the total magnetic flux 
up to 12^, the flux strip configuration turns out to be meta-stable. Likewise, fixing 
the value of the total magnetic flux at 12^, instability takes place when increasing the 
thickness up to /i = 7. As a result, two thick vortices of winding six form. Figure [T^ 
shows both configurations. Similar structures of intermediate states have been observed 
experimentally [IHl |16]. However the ones described in [15] and [16] do not arise as a 
consequence of a phase transition in the presence of an external magnetic field but as a 
result of the restoration of the symmetry in a type-I superconducting film as a uniform 
external auxiliary magnetic field H is increased. 

In order to test the meta-stability/instability of a flux strip, we proceed as follows. 
We first let the system relax to some relative minimum energy state compatible with 
periodic boundary conditions along the 2;-axis. That gives rise to a non-straight strip like 
the one in figure [TOT b). Both the electric field and the scalar field momentum are zero 
at that moment. Afterwards we implement small perturbations <^ |0| mm randomly 
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Figure 14. Superconducting film simulation, (a) Meta-stable profile of the scalar field 
along the junction at the meridional plane of the film y ^ 0. (b) Meta-stable profile 
of the magnetic fiux density along the strip at x — lx/'2- ^e.m. = 247r/e. For ^ = 3, 
twelve minima in \(f)\ and twelve maxima in By coexist, though no vortex forms. For 
/i = 7, two vortices of winding six form with the magnetic flux profile of equation (I63p . 



along the strip, once at a time. The size of the fluctuations in the transverse direction x 
vary from one lattice spacing up to the width of the strip. In the longitudinal direction 
z the size takes values from one to ttR. If the perturbation does not grow the system 
goes back to the original strip structure. That a strip is meta-stable and not stable can 
be easily seen by comparing its energy with the one in which the same amount of fiux 
is confined in separated stable vortices. The latter is energetically favourable. 

6. Formation of spontaneous topological defects 

In previous sections we have shown how local vortices and cosmic strings form in the 
collision of bubbles by trapping of an external magnetic fiux. A similar mechanism 
applies to the formation of topological defects when, instead of an external magnetic 
field, magnetic fiux is provided in the form of thermal fluctuations. One can think of 
the magnetic field long wavelength fluctuations as a uniform external magnetic field. 

The actual magnetic field fluctuations and nucleation of the actual bubbles need of 
a Monte Carlo simulation like the one performed in [18] . The analytical computation is a 
difficult task. In the finite-temperature field theory a first order phase transition effective 
potential of the form ([8]) comes from radiative corrections to the zero-temperature 
potential. In the particular case of the t/(l)-invariant renormalizable theory, 

£r=o = -^Ff^^F^" + D^4>D^'cf>* + A(#* - r/^)^ rj, X constants. (64) 
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The first order transition effective potential has been calculated in the type-I regime 
(A/e^ < 1) in several contexts. Some of them are concerned with quantum corrections 
[S] while others deal with thermal perturbations [IHl SHI HZl [IHl HE]. Because > A, 
radiative corrections due to the gauge field are dominant. The nucleation rate and the 
bubble structure are affected by the interaction of the fluctuation of both the scalar field 
and the gauge field [l9l [50] • In particular, this implies that the bounce solution and 
the nucleation rate as calculated in [301 EI] are only approximations. We will therefore 
simplify the problem by considering that the nucleation rate F, the bubble expansion 
rate v and the typical bubble nucleation radius Rn are free parameters which must be 
ultimately given by the actual field theory. 

Thus, the finite-temperature vacuum manifold consists of the set of minima of a 
potential like — A(00* —V^Y radiatively corrected. The equilibrium state corresponds to 
the Gibbs ensemble. Thermal fluctuations of the gauge field are strongly suppressed by 
Meissner effect in the broken phase where the photon is heavy. They are however present 
in the large areas of symmetric phase between bubbles. Therefore, with regards to the 
production of spontaneous magnetic field, we will concentrate on the long-wavelength 
fluctuations of Ai. 

6.1. Estimation of the typical winding number of thermal topological defects 

We will review here the arguments and predictions presented in [22] and [T8] . 

Let us assume that the dynamics of the whole system is adiabatic. Therefore, 
bubbles expand slowly enough such that the magnetic field can always be considered in 
thermal equilibrium in the symmetric phase areas between bubbles. Thus, the spectrum 
of the wavelength fluctuations corresponds to the long wavelength spectrum obtained 
in the thermal ensemble where the saddle point approximation is accurate enough given 
that the fluctuations are nearly gaussian. The spectrum is so the two-point correlation 
function of the gauge field in the symmetric phase. We have already calculated the 
functional form of that spectrum in the context of the zero-temperature field theory 
in section [51 In the thermal ensemble it is the Boltzmann factor exp {—F/T) what 
determines the weight of a configuration. Therefore, the spectrum of fluctuations is 
the Green's function of the equation of Ai for an extreme of the function F/T, that 
is: V^Gjj = TS'^{r — r')5ij. [j The label n denotes the number of space dimensions 
and so the functional form of the magnetic field fluctuations depends on the spatial 
dimensionality: 

(i?(A:)5(g)r ^T(27r)V(A:-g), (65) 
{By{k) By{q}f^' ^ T{27rmk - q}^-\kl (66) 

* This is a good approximation as long as only long wave length fluctuations are considered. It 
corresponds to the classical Ray leigh- Jeans spectrum of balckbody radiation and so valid for |A:| < T. 
In the case in hands we are considering l/|fc| > 2/m/i > 2/T. For the Lagrangian ([64l) . the strength of 
the last inequality depends on the coupling A and the v.e.v. of the scalar field. 
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{B,{k) B,iq))^^ ^ T{27if5\k + - (67) 

Formulas (165|) . (l66l) and (!67|) correspond to the fully two-dimensional, superconducting 
and fully three-dimensional models respectively. The superconducting one (|66ll was 
derived in following similar steps to those we followed in the computation of (15^ 
in the context of the zero-temperature theory. Again, top-arrows label two-dimensional 
vectors living in a superconductor film which is embedded in a three dimensional space 
and By stands for the component of the magnetic field orthogonal to the film. Generally, 
we will write the formulas above as 

{B,{k)BM)r = T(27r)"5"(A; + q)G^{k), (68) 

where the label m denotes the model to which the function corresponds and the indices 
i, j and the integer n take their values accordingly. 

As bubbles expand the magnetic flux associated to those fluctuations is expelled by 
Meissner effect. When several bubbles meet, a simply connected region of symmetric 
phase arises at the center of the collision. The magnetic flux in that region gets trapped 
by the mechanism explained in section 13.21 and a thick vortex or heavy cosmic string 
form. The typical amount of flux trapped can be estimated under adiabatic conditions. 
As the free energy F"^ is quadratic in the magnetic field components Bi{r), the magnetic 
energy F^^^/T = ^ J Bi{r)B^{r)d^r can be written in thermal equilibrium as 

Kag = ^/^"^/ d-r'B\r){G^)-\ry)B^{r'), (69) 

where (G'™')~-^(r, r') is the inverse of G^{k) defined in ( 168|) and so the two-point 
interaction potential. In the fully two-dimensional model the thermal fluctuations are 
not correlated, (G^'^)^^(r, r') = S'^{r — r'). In that case, the equilibrium value of (l69ll 
can be exactly computed on a surface of size 5* in function of the total magnetic flux 
M ■ That is 

<J)2 

Tp2D,S ^ ^E.M. /yQ\ 
mag 2,S ' V / 

In the superconducting model (SC) and the fully-three dimensional model, because the 
thermal fluctuations are correlated, the integral in fl69|) does depend on the shape of the 
surface S. Nevertheless the magnetic energy is a quadratic function of the total flux in 
both cases: 

f-r/ « (71) 

The factors of proportionality which are missing are geometrical coefficients of order 
one. 

As the free energy is quadratic in the magnetic flux in thermal equilibrium, the 



Topological Defects from First Order Gauge Theory Phase Transitions 35 



$2 



thermal fluctuations of ^e.m. adjust to the gaussian distribution ~ exp ( — in two- 
dimensions and to ~ exp (— ^^'^ ) in three dimensions. Correspondingly, the typical 
value of the thermal flux in a surface of area S is the root-mean-square of ^e.m.'- 

^T^l"" = (73) 
« VtS^/\ (74) 

The typical size S" of a region of symmetric phase enclosed by bubbles can be estimated 
using dimensional arguments. The characteristic length of that region must be of the 
order of the bubble radius. If the transition is strong first order, the radius is determined 
by the nucleation rate V and the expansion rate f . S* ~ (i;/r)^/^ holds in the fully 
two-dimensional and the superconducting models while 5* ~ (f/F)^/^ does in the fully 
three-dimensional model. 

Making use of the estimations above we can compute the typical winding number 
of vortices and cosmic strings in each case: 

Nl%,. oc f rV2(,;/r)V3, (75) 
« ^TV^(./r)V^ (76) 

NSy,. oc f TV2(./r)^/«. (77) 

A variety of windings is expected in the actual distribution. The properties of the 
winding numbers distribution function are to be reported in [3^ 



6.2. Flux trapping vs. Kibble mechanism 

As it was mentioned in section 13.21 (see also section 17.31) . Kibble mechanism can be 
used in the form explained there to explain the formation of local topological defects 
in absence of magnetic field. On the other hand, because the typical winding predicted 
by Kibble mechanism is of the order of the square root of the number of bubbles 
involved in a typical collision, it must be of order one. In particular, if the width 
of the bubble wall is negligible, it is exactly one. That is so because the probability for 
more than three bubbles to collide at a single point is zero. Therefore the trapping of 
thermal flux will be dominant over Kibble mechanism in the formation of vortices as 
,typ. > 1. In order to estimate when that is so we have to find out a relation between 
the parameters involved in the equations (1751) . fl76l) and fl77|) . First, we will assume that 
the electromagnetic coupling e is of order one. That is a reasonable assumption in both 
GUTs and electromagnetismj^l Then, we have to compare the value of the temperature 
T to the value of some power of ^ which is equivalent to some power of the radius of 
the bubbles. As mentioned at the beginning of this section, F is difficult to estimate 



f It is however still missing in this calculation the effect of the magnetic permeability in the 
superconducting model. 
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and so is f . We will take a more conservative approach. If the transition is not weak 
r is exponentially suppressed either by exp [—F{(f)T)/T] if the transition is temperature 
driven or by exp [—Se{4>o)] if it is a quantum transition. Se{4>o) stands for the Euclidean 
action of the bounce solution at zero-temperature while F{(j)T) does for the free energy 
of the bounce solution at finite-temperature T. Also, the size of the bubbles when they 
meet must be larger than or equal to the initial nucleation radius at the transition 
temperature R^iT). Following Linde's arguments [3T] . 

r ~ exp [-F(0t)/T] if F{<f)T)/T<SE{h), (78) 
r ~ exp [-Se{(Po)] if Se{(I)o) < F{<Pt)IT, i?jv(0) < {2T)-\ (79) 

If the transition is temperature driven, > T~^ . Therefore, the above conditions 
together with our estimates in equations ( 1751) . ( 1761) and ( 177|) imply that A''^; > 1 is 
guaranteed if the transition is strong first order and temperature-driven. For the theory 
in (!64l) that is the case in the couplings range > A > e^. Otherwise not only the 
nucleation process but also the bubble expansion rate must be considered. 



7. Josephson-Iike currents at the junctions and the geodesic rule in gauge 
theories 

7.1. The gauge-invariant phase difference 

In section 12.21 it was argued that it is D^j^O the field to be tracked in the gauge theory in 
order to describe the formation of topological defects. 

The definition of D^O in ffTOj) can be inferred from the application of the differential 
operator Z)^ onto the product exp {iO) ■ |0|. Following the Leibnitz's rule, 

D^{ei,.Y>{id) ■ 101) = exp(z^^)D^|0| +zexp(z^)|0|D^^. (80) 

On the other hand, the LHS of ( !80l) can also be written as 

D^(exp {lO) ■ 101) = {d^ + leA^M exp {lO). (81) 

By expanding the RHS of flSTl) and using the fact that the covariant derivative of a 

neutral function like |0| is equivalent to its ordinary derivative in (IHOl) . identity ( ITOj) 
holds. 

The local observables can be written in terms of the gauge- invariant quantities |0| 
and D^j^O as: 

Cooper pairs density |0p, (82) 

magnetic field B= (1/e) V x DO, (83) 

electric field E= (l/e) {doDO -VDoO), (84) 

charge density qe = 2e|0|^Z?o^5 (85) 

electrostatic potential Vei. = (1/e) DqO, (86) 

Noether current density j = — 2e|0|^D6', (87) 
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where in this case arrows denote spatial components vectors. 

Likewise, as d^dyO = dyd^O, it is possible to write down the equations of motion in 
terms of |0| and D^9 alone: 

[dl + adt - - D^OD.e + = 0, (88) 

d^d^D^e - d^'d^DJ + 2e^\(j)\'^D^e + adtD.e = 0, (89) 

where in the last equation we can identify ed^Fy^ = d^dyD^O — d^d^DyO and 
-eU^^+]f"') = '2e'\(p\'D^e + adtD.e. 

In the global theory D^O reduces to d^O and the phase difference dxO dx = 
ttB — ctA = Aaaf, is meaningful. In the gauge theory, when no magnetic field is present 
and the gauge = can be chosen, 'jab = /^^ D^O dx = Aaab as well. In this case, 
the only reason for the existence of non-zero D^j^O is the spatial overlap between two 
different vacuum domains A and B as it occurs in the junction between two colliding 
bubbles. In absence of magnetic field, D^9 can be interpreted as an induced spontaneous 
connection as a result of the coupling of the light ('fast') degrees of freedom living in 
the junction to the heavy ('slow') degrees of freedom living deep in the broken phase 
where = A/2e|0|mm- Even in the actual gauge theory this spontaneous connection 
cannot be gauged away. One can nevertheless perform the gauge transformation 

A, -> (1/e) dj, 

d^e (90) 

in such a way that DxO remains invariant and so does 7afe. 

If a magnetic field Bi does exist there is an extra divergenceless piece of connection 
induced by B^. It is always possible to write D^O as a sum of an exact derivative diO 
plus a divergenceless component eAi in the Coulomb gauge: 

DiO = diO + eAi such that d'Ai = 0. (91) 

As dD6 = edA = eF, did is the locally flat component of the connection while Ai carries 
the curvature in the gauge (19T!) . 

In summary, the local connection Di6 is composed of two parts: 
•A spontaneous, flat component did which, in the Coulomb gauge, carries the 
longitudinal modes of Di6. Its origin is in the breakdown of the U{1) symmetry. In the 
global theory, its interaction with the neutral scalar particles gives rise to the geodesic 
rule. 

•An induced, curved component eAi which, in the Coulomb gauge, carries the transverse 
modes of Di6. Its origin is in a net magnetic field. In the gauge theory, its interaction 
with the charged scalar particles gives rise to the Aharonov-Bohm effect. 

Both the local quantities diO and eAi and their integrals J did dx^ and e j Ai dx^ 
come to be respectively equal (up to a sign) deep in the broken phase where Di9 = 0. 
As one approaches the center of a vortex of winding at r = 0, Di9 is not zero any 
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more for r < ^ V^w^l, where the magnetic flux is confined. The holonomy of any 
loop C of radius r enclosing the center of the vortex, §(jdi9 dx"^ = 2ttN^, does however 
remain unaltered as r goes to zero by continuity. It is the discrepancy between the 
holonomy and the magnetic flux in the circle of radius what gives rise to the closed 
steady screening currents in the core of the vortex. 



7.2. DC Josephson-like effect across bubbles junctions and flux strips 

To our knowledge Weinberg was the first to show the periodicity of the Josephson 
currents using field theory arguments [51]. This will be our starting point to study the 
electric current traversing the junction between two bubbles as they collide. 

We will consider stationary situations and so neglect time oscillations, i.e. Dq9 ^ 0, 
dt\(f)\ ~ 0. Gauge invariance implies that, assuming no gradients along the junction 
between two superconducting pieces separated by a gap of width \b — a\, the value 
of the current self-interaction energy through the junction is an even function of 
'Jab = fa Dx9 dx: 

oo 

£ab =^Cm COS {m'Jab), (92) 
m=0 

where {cm} are constants which incorporate the degrees of freedom of |0| in a < x < 6. 
If stationary conditions hold, the current density is [5Tj 

If the current flow is parametrised with the coordinate x, a uniform steady current 
turns up in a < X < 6 ^33j if linear superposition of the scalar flelds is valid, that is, 
if 2/ m/j <^ 1 6 — a I . Uniformity of along the integration path a < x < b implies that 
the main contribution to formulas (]92|93p comes from the terms m = and m = 1. 
Uniformity of the reduced charge- carriers density |0|g implies that cq ~ ci ~ — 2|0|q/x. 
Therefore, under these conditions, fl92l) and fl93l) read respectiveljjj] 

eab = ^^[cos{jab)-l], (94) 
X 

— 2e 

Jaf, = |0|osin(7,6), (95) 

X 

where \(f)\o is given in fl28|) . The potential flMl) is a Josephson coupling and fl95l) is a 
direct Josephson-like current though its origin is not in quantum-tunneling [52]. Its 
amplitude, 

Jc^'^\cf>\l (96) 

is instead a function of the reduced charge- carriers density and the length scale x- 

In the actual collision of two superconducting bubbles or blocks, steady and 

I A more precise computation of this equation can be found in 33J. 



Topological Defects from First Order Gauge Theory Phase Transitions 



39 




Figure 15. Fully two-dimensional model simulation, (a) Transverse current density jx 
along the junction at the middle of the gap between two superconducting blocks. The 
setup corresponds to that of figure[8]with ab^aa — — 37r/2. Time t is set to zero at the 
moment the blocks get in touch. The current amplitude grows with time as the blocks 
overlap and |0|s=o increases. The current profile adjust to that in (|95p . (b) Current 
self- interaction energy £ab as a function of jab (equation ([M])). Sab is normalized to 



uniformity conditions approximately hold at the middle of the junction. The charge 
carriers density |0p is roughly uniform at the very middle of the gap, x = 0, where 
101 presents a minimum and so S^^l^l^-^o = (figure [T71(b)). The non-zero value of 
Dx6 concentrates also in a very narrow strip along the junction with center at x = 
and width When dz\(j)\ is negligible the current flows across the junction in the x- 
direction. Figure [T5l(a) shows the shape of the transverse current density j^^ in an actual 
simulation. It is in good agreement with expression ( |95l) . The current is measured at the 
middle of the gap {x = 0) between the two superconducting blocks of figure [HI There is 
a magnetic flux of value Gvr/e uniformly distributed in the gap so that the gradient of 7^6 
is uniform in the direction z along the junction. Figure [TST b) shows the corresponding 
current self-interaction energy as a function of ■jab normalized to -\4>\l = 1. 

An interesting scenario in which an actual steady superconducting current can 
be found is a meta-stable flux strip. In general, neither the density of charge carriers 
nor the magnetic flux are uniform in a meta-stable strip like the one in figure [TH with 
/i = 3. The gradients of |0| and D^O across the strip make equations (I94II95P inapplicable. 
Instead, the general expression (p2|) must be used. In addition, the gradients along the 
strip must be incorporated in the coefficients of the sum. Therefore, eab{z) takes the 

§ In ((27| . X is what the authors of [2^ identify with the width of the 'wave packet' D^O. In 33] it is 
shown that x is of the order of m^^- 
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general form 

oo 

£abiz) = ^ Cmiz) COS [ni'jabiz)]. (97) 



m=0 

Correspondingly, the transverse current density is 



e m Cmjz) sin [m'jab{z)]. (98) 



m=l 



The functions {Cm{z)} are not completely arbitrary. All the physical quantities are 
periodic along a meta-stable strip. If the length of the strip is Lz and it contains 
approximately winding units of magnetic flux, their period is Az ~ L^/N^. 
Consequently, both the functions {Cm{z)} and jabiz) have period Az along the junction 
(modulo 27r in the case of jab)- In terms of the total magnetic flux $, Az = 

A first consequence of the current density distribution fl98l) along a flux strip is 
an interference pattern of the transverse total current different to that in an actual 
Josephson junction. Taking into account the symmetry and periodicity properties of 
the functions in jab{z), it can be formally integrated along the junction of a meta-stable 
strip in a superconductor film of thickness /i and total flux $. That is 

Iab= / f^Jabiz)dz 

^, , , , rsin (me$/27r)i , , 

C:.sin(mAa.,)[— i--^J, (99) 

m=l ' 

where the coefficients C'^ are constants and Aaab stands for the spontaneous phase 
difference. The thickness fi is taken large enough such that j^^ is approximately uniform 
in the direction normal to the film. Figure [16] shows the functional form of the first 
three harmonics. 

The met a-st ability of a fiux strip requires the gap width |6 — a| to be relatively 
broad. Broad means in the first place that |& — a| is greater than l/mt so that the edges 
of the strip are approximately straight and the linear magnetic fiux density is nearly 
uniform. Broad also implies that / dx:\(t>\ dx <^ / 101 dx, x 1^ "~ '^\- If these 
conditions are satisfied the density of charge carriers is nearly uniform in the range of 
interest. Therefore the first harmonics m = 0, m = 1 are the dominant terms in the 
sums (EMn]). 

Because the total current only vanishes for integral values of the total fiux in 
units of 27r/e, there must be a net current passing through the gap during the process 
of formation of a meta-stable fiux strip. That total current stops when the fiux in the 
junction gets quantized. It is worth noting that the remaining current density jab{z) 
along the strip does not radiate because of topological reasons. It is the fiux quantization 
rather than the existence of a mass gap what prevents this current from dissipating. 

The agreement between the lattice simulation and the Josephson-like current in 
(!95|) deserves some comments. At the moment the blocks get in contact /5^/2 '^•^ 
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Figure 16. Interference pattern of the current across a flux strip. denotes the m*'' 
term of the series (jHH)- 



is of order vr and, as we mentioned in section HJ the dynamics of jab is governed by a 
sine-Gordon equation. In the numerical simulation, this apparent problem is overcome 
naturally with the use of the non-compact U{1) lattice formulation [53]. As described in 
Appendix B it is the use of the link field Uj = exp {ie6xAi) that incorporates the non- 
linear effects. Hence, the Josephson-like coupling —■^J2i[^^^*(x)^hi^)^iy:+i) ~^ I'/'Cx)!^] 
in (IB.Sp and the sinusoidal functions of figure [T^a). 



7.3. The geodesic rule in gauge theories 



In [M] it was argued that the geodesic rule in a gauge theory holds by equilibrium 
arguments, i.e. as a result of the minimization of the energy. In [5S] it was shown 
that if no magnetic field is present the geodesic rule holds in the gauge theory by 
dynamical arguments, result of the equations of motion. A more detailed study 

of the bubble collision dynamics was carried out in [26] in terms of gauge-invariant 
quantities. We will show with the aid of numerical simulations that both energetic 
[5l] and dynamical [26] arguments hold and agree with the geodesic-magnetic rule 
presented in section 13.21 The global geodesic rule as proved in [55] is only satisfied 
if some particular prescription for the gauge choice can be given (see below). 

Let us consider two superconducting bubbles A and B facing each other. Once 
more, if the radius of the bubbles is much larger than the bubble-wall thickness the 
contact area can be well approximated by two superconducting blocks with a gap in 
between. Block A is nucleated at the point a a while block B is at the point of 
the vacuum manifold. The gap in between is initially in the symmetric phase and no 
magnetic flux at all is in it. Let us take as an example Aoab = O-b — O-a = Sir/ 4 and 
let the two blocks expand and coalesce. Figure [TWa) shows the time evolution of both 
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Figure 17. Fully two-dimensional model simulation, (a) Time evolution of the 
gauge-invariant phase difference jab through the junction a < x < b between two 
superconducting blocks with spontaneous phase difference Aoab = 37r/4. (b) Time 
evolution of the scalar field profile \4>\ across the junction a < x < b. 



'^ab and A6ab across the junction a < x < b. In agreement with the geodesic-magnetic 
rule, 7afe decreases in time until vanishing. In the contrary, A6ab remains constant and 
equal to Sn/i. The complex phase 9 does interpolate between and 37r/4 as stated 
by the geodesic rule because the gauge = has been chosen initially. However, at 
the moment the blocks meet we could have made the following gauge transformation 
everywhere: 

A, ^ (1/e) dj, 

dj 0. (100) 

The complex phase would have been uniform in both the blocks and the junction. In 
that case the evolution of '-jah would have still been the one shown in figure [TTT a) while 
AOah would have remained equal to zero. This shows that the physical modes are those 
of D^O, neither the Goldstone modes dxO nor the longitudinal modes of the gauge field 
Ax- When the temporal gauge Aq is chosen, the evolution of d^O and A^ just depends 
on the initial fixing of the remaining gauge degrees of freedom. Thus, the geodesic rule 
as applied in the global theory is a gauge artifact in the gauge theory. It holds in the 
gauge theory when the gauge choice = can be made [55], which is only possible 
in the absence of a magnetic field. One can however consider that in the absence of a 
magnetic field topological defects form by Kibble mechanism. After all, in the gauge 
theory it is also the spontaneous breakdown of the symmetry that gives rise to the DiO 
modes. 

In energetic terms what is happening is that, as bubbles coalesce, the current tends 
to vanish as 7afc goes to zero. The current could also vanish if evolved towards 
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TT. However vr corresponds to a maximum and so unstable point of the current self- 
interaction energy Eat (figure [TST b)). The current self- interaction potential Eat does 
decrease as jab evolves toward its minimum at 7afe = 0. 

In dynamical terms, the Maxwell equation (!89|) becomes a generalized Klein-Gordon 
equation for D^9 in the gauge 

d^D^e = const., (101) 

[d^d^ + adt + 2e^\(P\^]D^e = 0. (102) 

Integrating the resultant equation for D^O across the junction (using the rough 
approximation |0| ~ uniform), the equation for is 

[dtt + adt + 2e^\(P\^]jab = 0. (103) 
The study of the time-evolution of 'jat was carried out in detail in [26] taking 

a(B) - a(A) = 3it/4 

jt/e 
37i/4e 

0) 

_2 

CO 

> ./2e 
_2 
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Figure 18. Fully two-dimensional model simulation. Time evolution of the magnetic 
flux associated to the fluxon ring around the collision axis of two bubbles. The 
spontaneous phase difference between the bubbles is Aaab = 37r/4. 

101 ~ |0|mm, dx\(t>\ ~ 0. However, because the typical phase difference between any 
two bubbles in a multi-bubble collision is of the order of vr for a vortex to form, (11031) 
is not a good approximation. It is instead the sine-Gordon equation (l27|l what should 
be used in the first stages of the collision. It is only when jab relaxes to jab ^ tt that 
(11031) is valid and jab goes roughly as ~ (a^ — a^)e The numerical simulation 

of figure fTWa) so confirms iil 

Figure [TS] shows the evolution of the magnetic fluxon ring generated around the 

I In the simulation of figure [TTf a) the dynamics of jab is overdamped so that oscillations in time are 
suppressed. 
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collision axis when two bubbles meet with — = 37r/4. Because the collision there 
is softjjj] the fluxon can spread out in the directions orthogonal to the collision axis. 
That is the reason for its constant asymptotic value. 

8. Conclusions and outlook 

In the Abelian Higgs model we have shown that it is the local dynamics of the gauge- 
invariant field Di6 = diO + eAi and the global dynamics of the gauge-invariant phase 
difference 7 = / DiO dx^ that determine the formation of vortices in superconductors 
and hypothetical cosmic strings in grand unified theories. 

In the mean-field theory approximation we have found that, when a first order 
phase transition takes place in an external magnetic field, a variety of objects of different 
winding form. These are meta-stable strips and membranes of magnetic flux, chains of 
single-winding vortices, thick vortices and heavy strings. Thick vortices turn out to be 
stable in thin superconducting films up to some critical winding value. The cause of 
the stability and/or meta-stability of vortices and strips is a Coulomb-like long-range 
magnetic interaction induced by the gauge field which propagates out of the film. It 
might determine the order of the phase transition of a thin superconducting film. A 
Monte Carlo simulation similar to the one in [18] with an external magnetic field is 
needed to investigate this matter. 

The spontaneous breakdown of the U{1) gauge symmetry manifests in the location 
of the center of mass of single-winding vortex chains and finite size Abrikosov lattices. 

Under steady conditions we have derived a direct Josephson-like current traversing 
the meta-stable flux strips. For that calculation we have invoked non-local effects within 
a length scale x which we identify with the correlation length of DiO. 

Thermal long-wavelength fluctuations of the gauge fields get confined in a variety 
of winding topological defects in a first order phase transition. In the particular case of 
the fully three-dimensional Abelian Higgs model, a multi-tension cosmic string network 
arise [18]. The typical winding of the strings is ~ ^T^/^(f /F)^/^ if the transition is 
strong first order, > A > e^. For a thin superconducting film the prediction is 
~ ^T^/^(f /F)^/^. Even though the latter is not easy to test experimentally, it is our 
belief that trapping of an external magnetic flux in the form of strips, vortex chains 
and thick vortices can be achieved experimentally by inducing nucleation of bubbles in 
a type-I superconducting film sample. 

In the context of brane-world cosmology and extra dimensions phenomenology, it 
would be interesting to study the hypothetical non-Abelian currents of grand unified 
theories traversing a 3 + 1 dimensional brane (with or without fluxes) in between two 
bulk disconnected - or weakly overlapped - vacuum domains. 

II If the collision is violent, the highly non-equilibrium dynamics of the scalar field can lead to the 
formation of unit-winding vortex-anti vortex pairs ^56! giving rise to an apparent violation of the 
geodesic rule. 
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Appendix A. Breakdown of the translational invariance of the setup in 
figure [8] 

Let us consider the fully two-dimensional theory and take two parallel superconduct- 
ing blocks separated by a finite gap as in figure [H The blocks are semi-infinite in the 
transverse direction x and infinite along z. A non-zero magnetic field exists in the gap, 
X G [—d/2,d/2], and it is invariant along the 2;-axis. Before the blocks get in contact, 
neglecting the tension of the walls, the initial equilibrium configuration for the magnetic 
flux corresponds to that of a uniform magnetic field B. If either the walls tension were 
considered or we were working in the framework of the superconducting model, transla- 
tion invariance along the junction would also hold but the magnetic field strength would 
show a gradient across the gap. 

We aim to show that, regardless of the choice of gauge: a) local translation invariance 
along the junction does hold as long as a finite width gap in which |0| = V^; exists 
between the two blocks, b) local translation invariance breaks down as the superconduct- 
ing blocks get in contact. Condition a) is necessary for both SSB and for the stability 
of the straight strip in figure [9t^a). Condition b) is necessary for the formation of the 
structures presented in figures [9t^b) and[9](c). 

Proof of part a). 

The initial conditions of the local physical quantities are: 

1) Translation-invariance along the z direction: 

dMt=o = o Vx, 

d,B\t=o = Vx. 

The above equations also imply: 

dzf \t=o =0, at X such that dz\(p\x,t=o = 0, 
dzD'^6\t=Q = 0, at X such that dz\(t)\x,t=Q = and |0| > 

2) The system is initially at rest: 

dt(t)\t=o = Vx, 
Ei\t=o =0 Vx, 
dtj%=o = Vx. 



(A.l) 
(A.2) 

(A.3) 
(A.4) 
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3) Asymptotic vacuum state conditions: 

\4>\t=o = \4>\min for X < -d/2 or X > d/2, 

Di9\t=o = 0, for X < -d/2 or x > d/2. (A.5) 

4) Existence of a narrow gap of symmetric phase at the center: x G (— e, e) Wz, < e <^ d. 

|0|t=o = 0, for -e < X < e. (A. 6) 

Bearing in mind that all the fields are evaluated at t = 0, we will drop the subscript 
t = 0. Breakdown of translational invariance along z requires that 

d,{d^\^{x,z)\) ^ for some x G [~d/2,d/2]. (A.7) 

The equation of motion for |0| in the interval x G [—d/2, d/2] at t = subject to the 
above conditions is 

[d^ - + Die + Die + ^]|0| = 0. (A.8) 
Taking the derivative of (1A.8I1 along the direction z, 

d^dlicj^l - dldM - dM + W-^^ + d,{Dle + Diem = o. (a.9) 



The second, third and fourth terms in the LHS of equation ( lA.QIl vanish according to 
(lA.ip and we are left with 



d,dl\cf>\ = -d^iDle + DleM]. (A.io) 



In order to satisfy the inequality (lA.7p . the RHS of (lA.lOp must be different to zero at 



some point x G [—d/2, d/2]. The Maxwell equation (l89l) in the gauge of equation (110 II) 
reads 

[{dl + dl) - ml^{x)]D,e = 0, where m^^{x) = 2e^\(f){x)\\ (A.ll) 
Let us write the 0(2) massive field Die in terms of the polar coordinates g and ip: 

= Die + Die, 

DJ = gcos^ip), 

D,e = Qsm{i)), (A.12) 

so that the inequality (IA.7P implies 

dzdl\<P\ ^0 
^-9,(^^101)^0 

^ -2\(t)\Qd,Q^Q. (A.13) 

In ( jA.131) . Q does not vanish in the junction because there is a magnetic field. The value 
of |0(x)| equals in x G (— e,e) according to (lA.Gp . 2;-invariance of both |0| and B im- 
plies that dzQ = if 101 7^ (equations (lA.ll) . (]A.3I) ). Hence, the inequality (1A.13I) does 
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not hold and 2;-invariance remains as long as e > 0. Figures IXT] (a) and lAl( c) illustrate 
this situation. It is the case for a separation between blocks d ^ 2/m/j. 

Proof of part b ). 

Let us assume now that the blocks are not separated but they overlap instead in a 
region of width 2e <^ d. Thus, |0| > everywhere but, at most, at a finite number of 
points in the junction. Likewise the existence of a net magnetic field in the junction 
implies that g > in —d/2 < x < d/2. Therefore, in order to satisfy (]A.13I) we are left 
with the proof of dzQ 7^ at some x G [—d/2,d/2\. z-invariance holds at least in the 
range x < — e, x > e and therefore flA.ip reads now 



9^|0|j=o = X < -e, X > e. (A.14) 

We aim to show that, if a z-invariant magnetic field is present in the junction, z- 
invariance of |0| in the interval x G (— e, e) is not possible. 2;— invariance of the magnetic 
field can be expressed as a constraint equation for DiO, 

dlDJ = dAD,e, (A.15) 

and ( lA.llI) in the variables ( 1A.12I) reads 

[dl + dl - {dx^f - {d,i)f - mi{x)\Q = 0. (A.16) 



Let us write (1A.16) as 



[dl + dl-m'^'ix)]0 = O, (A.17) 

where we have defined ni^{x) = {d^i^Y + {dz'4^Y + '^7(^)5 and let us integrate (1A.17P 
from X —00 to X ^ +00. Because g falls deep inside the superconducting blocks as 
X goes to ±00, 

/+00 
dlgdx = dxg\+^ = 0. (A.18) 
-00 

Therefore, 

"+00 f+00 

52„ / '2 



d^gdx= I m^{x)gdx. (A.19) 

) J —00 

The LHS of ( 1A.19I) is greater than zero because the integrand of the RHS, m'^{x)g, is 
positive definite and both m'^{x) and g are greater than zero at some point in the gap 
for there is a net magnetic field and blocks overlap. That implies that, correspondingly, 
there must be at least some point xq in the interval [—d/2, d/2] at which d'^g\xo > 0. If 
this is so and dzg\xo 7^ as well, our proof is done and breakdown of the translational 
invariance of both \(f)\ and g is shown. If however dzg\xo = 0, then xq is a critical point. 
Because g is gauge-invariant it is continuous in the broken phase. Therefore there must 
be a neighborhood centered at Xq denoted by Q^o iii which dzg\n^^^ 7^ 0. Thus, the proof 
is complete. □ 

Figure lAll shows the breakdown of the 2;- invariance of |0| and = + j'f as the 
two blocks get in contact and there is a magnetic flux of value Gtt/c in the junction. 



Topological Defects from First Order Gauge Theory Phase Transitions 



48 




Figure Al. Fully two-dimensional model simulation, (a) t = 0. Initial profile of |(/)| 
according to the setup of figure Invariance along z holds, (b) t — 30. z-invariance 
breaks down as blocks meet (d ~ 2/m/i) and three minima (= b^e.m./'^t^) show up. 
(c) t = 0. Initial profile of the square of the current density = + j^. Invariance 
along z holds. The current density concentrates at the edges of the blocks where the 
magnetic field initially migrates due to small structure effects, (d) t — 30. z-invariance 
breaks down as blocks meet, gets concentrated around the minima of |0|. 



Appendix B. Lattice disretization 

To write down the discretized equations of motion, we define the forward and backward 
derivatives 

A^/W = ±5x-^ -/(.)) , (B.l) 

where 5x is the lattice spacing, z is a unit vector in the i direction and vectors are 
denoted by bold-faced letters. Similarly, we define the time derivative 

A*/(t) = 5t-' (/(i) - f^t-st)) , (B.2) 
where 5t is the time step. Using the link variable 



Uj = exp {ie6xAi) , 



(B.3) 
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we also define the corresponding covariant derivatives in the superconductor bulk 




(B.4) 



We make use of the Hamiltonian formalism to write the equations of motion ([6]), ([7]) as 
first order time-differential equations, '^'{t,^) is the conjugate momentum of the scalar 
field 0(t,x) and the electric field Ej_((_x) is the conjugate momentum of the gauge field 
Aj (f x)- In the temporal gauge Ao,(t,x) = 0, the equations of motion in the fully two- 
dimensional model are 

AtEi,(t,x) = ^ eijefczAjAfc A/^(t^x) — o'Ej_(t_x) 
jkl 

- 2eIm0^i_^)D+0(j,x), 
At7r(t,x) = Yl DrD,^0{t,x) - fT7r(t,x) 

i 

3 

- m]l(t>{t,x) + - fi;|0(t,x)|0{t,x) - 2A| 0(t,x)|V(t,x), (B.5) 

where ejj is the Levi-Civita tensor in 2D and the indices run from 1 to 2 in the 
summations. 

Analogously in 3D, 

AtAj^(t^x) = — Ej_(t_5t^x), 

jklm 

- 2eIm0^i_^)D+0(t,x), 
At7r(t,x) = D,"D+0(t,x) - fX7r(j,x) 

i 

3 

- m]l(t>{t,x) + - fi;|0(t,x)|0(t,x) - 2A| 0(t,x)|V(t,x), (B.6) 

where eijk is the Levi-Civita tensor in 3D and the indices run from 1 to 3 in the 
summations. 

In the superconducting model, the boundary conditions in f|T3l) are implemented 
modifying the equations for At7r(t^x) and A^E^ x) at the surface of the film. In the 
equation for /S.t'K^t,^) the index i in the summation only takes values 1 and 3 for y < — /u/2 
y > ji/2. Likewise, the term — 2elm0*^ ^^D+0(t x) is excluded from the equation of 
AjEy (t x) ioT y > ii/2 or y < — yu/2. The reason for these modifications is that the 
covariant derivatives as defined in the lattice, D+</)(x) and D~0(x), are ill-defined at 
y = /i/2, — /i/2 — 1 and y = — /i/2,/i/2 + 1 respectively. They must be discarded 
explicitly at those points because their vanishing is not possible in the superconducting 



Topological Defects from First Order Gauge Theory Phase Transitions 

phase according to formula (lB.4p . 

Once the above modifications have been implemented and the conditions 

0(t=o,x) = for y > fi/2 or y < -fi/2, 
a =0 ioT y > /i/2 or y < -/i/2 V t 

set up, the equations of motion in (1B.6I) take care of the conditions in (fT3l) . 
The lattice version of the Hamiltonian is, in 2D, 



tt2D 

^tot 2 



^*,{x)<P(x+5) 



-|0(x)P] 

1,3 = 1,2. 



In three dimensions. 



Hot 



E 



TT TT 



jk 



5x' 



5^[Re0^,)U, 



*,(x)<P(x+i) 



I0(x)| 
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(B.7) 



+ m^|0(,)|2 - + A|0(,)|^J , ^, J, = 1, 2, 3. (B.8) 

When measuring the energy associated to a superconductor film, the index i in the 
last summation does not take the value 2 associated to the y direction orthogonal to 
the film at ?/ = /i/2. It does however aX. y = — yu/2. 



Table CI. Figure number, coupling constants, lattice parameters and other setup values. 



fie;. no. 


e 


rriH 




A 


/ X / X / 


5x 


6t 


a 


^ini 


other parameters 


El 


0.2 


0.656 


0.295 


0.043 


384 X 384 


1.0 


0.01 


1.0 


0.0 




a 


0.2 


0.656 


0.295 


0.043 


384 X 384 


1.0 


0.01 


2.6 


6 • 10"^/e 




El 


0.34 


0.656 


0.295 


0.043 


384 X 384 X 384(p6c) 


1.0 


0.01 


0.25 


0.012/e 


d = 143 


H 


1.0 


0.656 


0.295 


0.043 


900 X 900 


1.0 


0.01 


0.0 


7.76 • 10"^/e 


d = 300 Rn = 14 




0.1 


0.656 


0.295 


0.043 


160 X 27r/? (, 


1.0 


0.01 


1.0 


3.0/e 
160xi?„,6,<, 


\b — a\ = 30 


h nl 
IIUI 


U.i 


U.DOD 


u.zyo 


U.U4o 


ZUU X oUypoc) X W()[poc) 


i.U 


U.Ui 


U.O 


n nocQ / ^ 
U.UZao/ e 


|o — a| = 4U, fl = o 


M 


0.1 


0.656 


0.295 


0.043 


384 X 384 


1.0 


0.01 


0.5 


0.011/e 


\b-a\= 30 




0.2 


0.656 


0.295 


0.043 


156 X 80 X 50{pbc) 


1.0 


0.01 


0.5 


0.1885/e 


ds = 10 


m 


0.2 


0.656 


0.295 


0.043 


156 X 80 X 50{pbc) 


1.0 


0.01 


0.5 


variable 


ds, IJL variable 


M 


0.2 


0.656 


0.295 


0.043 


160 X 50{pbc) X 50(p6c) 


1.0 


0.01 


0.5 


5.03- lO-^/e 


\b-a\= 30 




0.1 


0.2 


0.0151 


0.0014 


300 X 60{pbc) 


1.0 


0.01 


1.2 


0.0 


\b-a\= 30 


m 


0.2 


0.656 


0.295 


0.043 


384 X 384 


1.0 


0.01 


1.8 


0.0 


\b-a\= 30 


ED 


0.1 


0.656 


0.295 


0.043 


160 X 51(p6c) 


1.0 


0.01 


1.0 


0.0123/e 


\b-a\ = 30 



Appendix C. Table of parameter values in the numerical simulations 

,where x 1^ x ly is the lattice size in lattice spacing units; (p6c) stands for 'periodic boundary conditions'; 5x is the lattice spacing; 
5t is the time step; a is the damping rate; B^^i is the initial magnetic flux density; /i is the film thickness in the superconducting 
model; d stands for bubble centers separation; dg is the initial square side; is the initial bubble radius; \h — a\ is the gap size 
between blocks or bubbles. The value for a in fig.no. [7] is 0.0 in the equation of the gauge field and is variable in the equation for cj) 
as it determines the bubble expansion rate. 
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